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Abstract 

In this paper, we propose a kernel method for exact tail asymptotics of a random walk to neigh¬ 
borhoods in the quarter plane. This is a two-dimensional method, which does not require a deter¬ 
mination of the unknown generating function(s). Instead, in terms of the asymptotic analysis and a 
Tauberian-like theorem, we show that the information about the location of the dominant singularity 
or singularities and the detailed asymptotic property at a dominant singularity is sufficient for the 
exact tail asymptotic behaviour for the marginal distributions and also for joint probabilities along 
a coordinate direction. We provide all details, not only for a “typical” case, the case with a single 
dominant singularity for an unknown generating function, but also for all non-typical cases which 
have not been studied before. A total of four types of exact tail asymptotics are found for the typical 
case, which have been reported in the literature. We also show that on the circle of convergence, an 
unknown generating function could have two dominant singularities instead of one, which can lead 
to a new periodic phenomena. Examples are illustrated by using this kernel method. This paper can 
be considered as a systematic summary and extension of existing ideas, which also contains new and 
interesting research results. 

Keywords: random walks in the quarter plane; stationary distribution; generating function; kernel 
methods; singularity analysis; exact tail asymptotics; light tail 


1 Introduction 

Two-dimensional discrete random walks in the quarter plane are classical models, that could be either 
probabilistic or combinatorial. Studying these models is important and often fundamental for both the¬ 
oretical and applied purposes. For a stable probabilistic model, it is of significant interest to study its 
stationary probabilities. However, only for very limited special cases, a closed-form solution is available 
for the stationary probability distribution. This adds value to studying tail asymptotic properties in 
stationary probabilities, since performance bounds and approximations can often be developed from 
the tail asymptotic property. The focus of this paper is to characterize exact tail asymptotics. Specif¬ 
ically, we propose a kernel method to systematically study the exact tail behaviour for the stationary 
probability distribution of the random walk in the quarter plane. 

The kernel method proposed here is an extension of the classical kernel method, first introduced by 
Knuth [28] and later developed as the kernel method by Banderier et al. [3]. The standard kernel method 
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deals with the case of a functional equation of the fundamental form K{x,y)F{x,y) = A{x,y)G{x) + 
B{x,y), where F{x,y) and G{x) are unknown functions. The key idea in the kernel method is to find 
a branch y = yQ{x), such that, at {x,yo{x)), the kernel function is zero, or K{x,yo{x)) = 0. When 
analytically substituting this branch into the right hand side of the fundamental form, we then have 
G{x) = —B{x,yo{x))/A{x,yo{x)), and hence, 


F{x,y) 


-Ajx, y)B{x, yo{x))/A{x, yojx)) + B{x, y) 
K{x,y) 


However, applying the above idea to the fundamental form of a two-dimensional random walk does 
not immediately lead to a determination of the generating function P{x,y). Instead, it provides a 
relationship between two unknown generating functions 7ri{x) and 7r2{y), referred to as the generating 
functions for the boundary probabilities. This is the key challenge in the analysis of using a kernel 
method. Therefore, a good understanding on the interlace of these two functions is crucial. 

Following the early research by Malyshev [Ml SO], the algebraic method targeting on expressing 
the unknwon generating functions was further systematically updated in Fayolle, lasnogorodski and 
Malyshev mi based on the study of the kernel equation. The authors indicated in their book that: 
“Even if asymptotic problems were not mentioned in this book, they have many applications and are 
mostly interesting for higher dimensions.” The proposed kernel method in this paper is a continuation of 
the study in P. Research on tail asymptotics for various models following the method (determination 
of the unknown generating function(s) first) of [11] or other closely related methods can be found in 
Flatto and McKean CD, Fayolle and lasnogorodski [9], Fayolle, King and Mitrani m, Cohen and 
Boxma |7] , Flatto and Hahn [T5] , Flatto [H] , Wright [55] , Kurkova and Suhov [30] , Bousquet-Melou [6] , 
Morrison |50j . Li and Zhao |34l 135] . Guillemin and Leeuwaarden [20], and Li, Tavakoli and Zhao [32j . 

Different from the work mentioned above, which requires characterizing or expressing the unknown 
generating function, such as a closed-form solution or an integral expression through boundary value 
problems, the proposed kernel method only requires the information about the dominant singularities 
of the unknown function, including the location and detailed asymptotic property at the dominant 
singularities. Because of this, the method makes it possible to systematically deal with all random 
walks instead of a model based treatment. In a recent research, Li and Zhao [36] applied this method to 
a specific model, and Li, Tavakoli and Zhao [32] to the singular random walks. For exact tail asymptotics 
without a determination of the unknown generating function(s) or Laplace transformation function(s), 
different methods were used in the following studies: Abate and Whitt [I], Lieshout and Mandjes m, 
Miyazawa and Rolski [58], Dai and Miyazawa [8]. 

Other methods for studying two-dimensional problems, including exact tail asymptotics, also exist, 
for example, based on large deviations, on properties of the Markov additive process (including matrix- 
analytic methods), or on asymptotic properties of the Green functions. References include Borovkov 
and Mogul’skii [5], McDonald m, Foley and McDonald [m [n [IS], Khanchi [Ml [25], Adan, Foley 
and McDonald [2], Raschel [52], Miyazawa [441 [55l 146] . Kobayashi and Miyazawa [26], Takahashi, 
Fujimoto and Makimoto [53], Haque m, Miyazawa [43], Miyazawa and Zhao [49], Kroese, Scheinhardt 
and Taylor [29], Haque, Liu and Zhao [22] . Li and Zhao [33], Motyer and Taylor m, Li, Miyazawa 
and Zhao m, He, Li and Zhao [23], Liu, Miyazawa and Zhao [38] . Tang and Zhao [54], Kobayashi, 
Miyazawa and Zhao [23, among others. For more references, people may refer to a recent survey on 
tail asymptotics of multi-dimensional reflecting processes for queueing networks by Miyazawa m- 

The main focus of this paper is to propose a kernel method for exact tail asymptotics of random 
walks in the quarter plane following the ideal in based on which a complete description of the 
exact tail asymptotics for stationary probabilities of a non-singular genus 1 random walk is obtained. 
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We claim that the unknown generating function vri(x), or equivalently, 7r2(y), has either one or two 
dominant singularities. For the case of either one dominant singularity, or two dominant singularities 
with different asymptotic properties, a total of four types of exact tail asymptotics exists: (1) exact 
geometric decay; (2) a geometric decay multiplied by a factor of (3) a geometric decay multiplied 

by a factor of and (4) a geometric decay multiplied by a factor of n. These results are essentially 

not new (for examples see references O [T71 dU 051 [25]) except that the fourth type is missing from 
previous studies for the discrete random walk, but was reported for the continuous random walk in [8]. 
For the case of two dominant singularities with the same asymptotic property, a new periodic phenomena 
in the tail asymptotic property is discovered, which has not been reported in previous literature. For 
the tail asymptotic behaviour of the non-boundary joint probabilities along a coordinate direction, a 
new method based on recursive relationships of probability generating functions will be applied, which 
is an extension of the idea used in [36| . 

For an unknown generating function of probabilities, a Tauberian-like theorem is used as a bridge 
to link the asymptotic property of the function at its dominant singularities to the tail asymptotic 
property of its coefficient, or in our case, stationary probabilities. This theorem does not require the 
monitonicity in the probabilities, which is required by a standard Tauberian theorem and cannot be 
verified in general, or Heaviside operational calculus, which is usually very difficult to be rigorous. 
However, the price paid for applying the Tauberian-like theorem requires more in analyticity of the 
function and detailed information about all dominant singularities, or singularities on the circle of 
convergence. Therefore we need to provide information about how many singularities exist on the circle 
of convergence and their detailed properties, such as the nature of the singularity and the multiplicity 
in the case of the pole, for the random walk. It is not always true that only one singularity exists on 
the circle of convergence. Technical details are needed to address these issues. 

The kernel method immediately leads to exact tail asymptotics in the boundary probabilities, in 
both directions, based on which exact tail asymptotics in a marginal distribution will become clear. 
However, it does not directly lead to exact tail asymptotic properties for the joint probabilities along 
a coordinate direction, except for the boundary probabilities as mentioned above. Therefore, further 
efforts are required. In this paper, we propose a method, based on difference equations of the unknown 
generating functions, to do the asymptotic analysis, which successfully overcomes the hurdle for exact 
tail asymptotics for joint probabilities. 

The rest of the paper is organized into eight sections. In section 2, after the model description, 
the so-called fundamental form for the random walk in the quarter plane is provided, together with 
a stability condition. Section 3 contains necessary properties for the two branches (or an algebraic 
function) defined by the kernel equation and for the branch points of the branches. These properties 
are either directly from m or further refinements. Section 4 consists of six subsections for the purpose 
of characterizing the asymptotic properties of the unknown generating functions and TT 2 {y) at 

their dominant singularities. Specifically, two Tauberian-like theorems are introduced in subsection 1; 
the interlace between the two unknown generating functions is discussed in subsection 2, which plays a 
key role in the proposed kernel method; detailed properties for singularities of the unknown generating 
functions are obtained in subsections 3-5, which finally lead to the main theorem (Theorem 14.8p in this 
section provided in the last subsection. In Section 5, asymptotic analysis for the boundary generating 
functions is carried out, which directly leads to the tail asymptotics for the boundary probabilities 
in terms of the Tauberian-like theorem. In Section 6, based on the asymptotic results obtained for 
the generating function of boundary probabilities in the previous section, and the fundamental form, 
exact tail asymptotic properties for the two marginal distributions are provided. Exact tail asymptotic 
properties for joint probabilities along a coordinate direction is addressed in Section 7, which is not a 
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direct result from the kernel method. Instead, we propose a difference equation method to carry out 
an asymptotic analysis of a sequence of unknown generating functions. The last section contributes to 
some concluding remarks and two examples by applying the kernel method. 


2 Description of the Random Walk 

The random walk in the quarter plane used in this paper to demonstrate the kernel method is a reflected 
random walk or a Markov chain with the state space = {(m, n);m, n are non-negative integers }. 
To describe this process, we divide the whole quadrant into four regions: the interior = 
{(m,n);m,n = 1,2,...}, horizontal boundary = {{m,0);m = 1,2,...}, vertical boundary S 2 = 
{(0, n); n = 1,2,...}, and the origin Sq = {(0,0)}, or Z^ = 5+ U S*! U ^2 U Sq. In each of these regions, 
the transition is homogeneous. Specifically, let X_|_, Xi, X 2 and Xq be random variables having the 
distributions, respectively, pij with i,j = 0, ±1; with i = 0,±1 and j = 0,1; p\^j with i = 0,1 

and j = 0, ±1; and pf'j with i,j = 0,1. Then, the transition probabilities of the random walk (Markov 
chain) Lt = L 2 {t)) are given by 

PiLt+i = {m 2 ,n 2 )\Lt = (mi,ni)} = 

f P(X+= (m 2 - mi,re 2 - ni)), if (m 2 , 712 ) G 5, (mi, m) G 5+, 

\ P{Xj: = (m 2 — mi, 772 — ^ 1 ))) if im 2 ,n 2 ) G S, (mi, ni) G Sk with k = 0,1, 2. 

2.1 Ergodicity conditions 

A stability (ergodic) condition can be found in Theorem 3.3.1 of Fayolle, lasnogorodski and Maly¬ 
shev P, which has been amended by Kobayashi and Miyazawa as Lemma 2.1 in [26]. This condition 
is stated in terms of the drift vectors defined by 

M = {MxMy) = 

^ i j j i ^ 

M(‘) = {m<;>m»')=(e<ej>s)-5:^}Ej>s))’ 

^ i j j i ' 

mp) = {MfMf)=(5:i(5:pg>).5:i(5:pg>)). 

^ i j j i ^ 

Theorem 2.1 (Theorem 3.3.1 in |11| and Lemma 2.1 in |26| 1 When M / 0, the random walk 
is ergodic if and only if one of the following three conditions holds: 

1. <0, My <0, < 0 and MyM^ - < 0; 

2. Mx < 0, My > 0, MyM^'^ — MxMy"^^ < 0 and < 0 if M^'^ = 0; 

3. Mx > 0, My < 0, MxM^^ — MyMx^^ < 0 and M^'^ < 0 if M^^ = 0. 

Throughout the paper, we make the following assumption, unless otherwise specified: 

Assumption 1 The random walk Lt is irreducible, positive recurrent and aperiodic. 
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Under Assumption [H let 'Km,n be the unique stationary probability distribution of the random walk. 


Remark 2.1 It should be noted that for a stable random walk, the condition M ^ 0 is equivalent to 
that both sequences {vro,n} o,re light-tailed (for example, see Lemma 3.3 of which is not 

our focus of this paper. Therefore, Theorem \2.1\ provides a necessary and sufficient stability condition 
for the light-tailed case. 


2.2 Fundamental Form 

Define the following generating functions of the probability sequences for the interior states, horizontal 
boundary states and vertical boundary states, respectively, 

OO OO 

m=l n=l 

OO 

m=l 

OO 

71-2 (y) = 

n=l 


The so-called fundamental form of the random walk provides a functional equation relating the three 
unknown generating functions 7r(x,y), vri(x) and T^ 2 {y)- To state the fundamental form, we define 


h{x,y) 


hi{x,y) 


/i 2 (x,y) 


ho{x,y) 


/II 

\i=-ij=-i 

a{x)y‘^ + b{x)y + c{x) = d{y)x^ + b{y)x + c(y), 

x(y^ ^pQx*y^ - 

V *=-1 i =0 / 

ai(a;)y bi{x) = ai{y)x^ + bi{y)x + ci{y), 

\i=0j=-l ) 

a2{y)x + b2{y) = a2{x)y‘^ + b2{x)y + C2(x), 

\ i =0 j =0 / 

ao{x)y + bo{x) = ao{y)x + bo{y), 
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where 


o(x) = p_i,i+po,ia: 
b{x) = p-ifi-{l-pofi)x+pifix'^, 
c{x) = p-i-i+po-ix+pi-ix"^, 
a{y) = pi-i+pi,oy+pi,iy‘^, 

Hy) = Po-i - (1 - po,o)y + po,iy^, 
c{y) = p-i-i+p-i,oy + p-i,iy‘^, 

ai{x) = 1 + pga: + bi{x) = ^ “ (l “ Pofl) ^ 

ai{y) = +pS^Jy,6i(y) - 1+pSy,ci(y) =p5 o+P-J,iy 

02 ( 3 ;) = Po} + pf}x, b 2 {x) = p^Q^Q - 1 + pf^^X, C 2 {x) = 

a2{y) = pfh + pf^y + P?}y‘^,b2{y) = p^oh “ “ Pcho) v + 

ao{x) = pI)°} + pS°Jx, bo(x) = pf^^x - (l - Po^) , 

My) = pffi + p?}y,My) = Po}y - • 

The basic equation of the generating function of the joint distribution, or the fundamental form of 
the random walk, is given by 

- h{x, y)7r{x, y) = hi{x, y)'Ki{x) + h 2 {x, y)My) + Mx, y)'Xo,o. (2.1) 

The reason for the above functional equation to be called fundamental is largely due to the fact that 
through analysis of this equation, the unknown generating functions can be determined or expressed, 
for example, through algebraic methods and boundary value problems as illustrated in Fayolle, lasno- 
gorodski and Malyshev [H]- The kernel method presented here also starts with the fundamental form, 
but without expressing generating functions first. 

Remark 2.2 The generating function 7r(x, y) is defined for m,n > 0, excluding the boundary probabil¬ 
ities. h2.1\) was proved in (1.3.6) in fll]/ . Based on 112.1\} . one can also obtain a similar fundamental 
form using generating functions including boundary probabilities: n(x, y) = X]m=o ^’n=o '^m,nX^y^, 
ni(a;) = Em=o and U 2 {y) = E^o^o,ny’^- 

For the conclusion of this section, we can easily check the following expressions, some of which will 
be needed in later sections: 

My = a(l) - c(l) = S'(l) + 6'(1) + d'(l), = a{l) - d(l) = a'(l) + b'{l) + c'(l), (2.2) 

M(i) = ai(l) =h;(l) + 6;(l) +2^(1), mW =Si(l) -Fi(l) = a;(l) + 6;(1), (2.3) 

Mf = a 2 (l) - C 2 (l) = h'(1) + 6^(1), Mi^) = S 2 (l) = o'(1) + 6'2(1) + 4(1). (2.4) 

3 Branch Points And Functions Defined by the Kernel Equation 

The property of the random walk relies on the property of the kernel function h and functions hi and 
/i 2 . The kernel function plays a key role in the kernel method. 
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Definition 3.1 A random variable is called non-singular if the kernel funetion h{x,y), as a polynomial 
in the two variables x and y, is irreducible (equivalently, if h = fg then either f or g is a constant) and 
quadratic in both variables. 

Throughout the paper unless otherwise specified, we make the second assumption below. 

Assumption 2 The random considered is non-singular. 

The non-singular condition for a random walk is closely related to the irreducibility of the marginal 
processes Li[t) and L 2 {t), but they are not the same concept. A necessary and sufficient condition for a 
random walk to be singular is given, in terms of pij, in Lemma 2.3.2 in m- Study on tail asymptotics 
for a singular random walk is either easier or similar to the non-singular case, which can be found in 
Li, Tavakoli and Zhao [32]. 

The starting point of our analysis is the set of all pairs {x, y) satisfy the kernel equation, or 

B = {{x,y) e : h{x,y) = 0}, 

where C is the set of all complex numbers. The kernel function can be considered as a quadratic form in 
either x or y with the coefficients being functions of y or x, respectively. Therefore, the kernel equation 
can be written as 

a{x)y'^ -h b{x)y + c{x) = a{y)x^ + h{y)x + c{y) = 0. (3.1) 

For a fixed x, the two solutions to the kernel equation as a quadratic form in y are given by 

_ -b{x)±y^Diix) 

’ 2a{x) 

if a{x) / 0, where Di{x) = h‘^{x) — Aa{x)c{x). Notice that non-singularity implies that a[x) ^ 0 and, 
therefore, only up to two values of x could lead to o(x) = 0 since o(x) is a polynomial of degree up to 2. 

Similarly, for a fixed y, the two solutions to the kernel equation as a quadratic form in x are given 
by 

-Kv) VB)2{y) 

2a{y) 

where D 2 {y) = b‘^{y) - 4a(y)c(y). 

It is important to study the set B, or equivalently Y±{x) or X±{y), since for all {x,y) G B with 
\TT(x,y)\ < oo, the right hand side of the fundamental form is also zero, which provides a relationship 
between the two unknown generating functions tti and ^ 2 . In the above, y^Di(x) is well-defined if 
Di{x) > 0 and similarly D 2 [y) is well-defined if D 2 {y) >0. As a function of a complex variable, the 
square root is a two-valued function. To specify a branch, when z is complex, ^/z is defined such that 
^=1. 

Let z = Di{x). Then, both y_(x) and b + (x) are analytic as long as z ^ (—oo,0] and a{x) / 0. 
For these two functions, we start from a region, in which they are analytic, and consider an analytic 
continuation of these two functions. In this consideration, the key is the continuation of y^Di(x). 

Definition 3.2 A branch point ofY±{x) (X±(y)J is a value of x (y) such that Di{x) = 0 (D 2 {y) = Oj. 
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To discuss the branch points, notice that the discriminant Di {D 2 ) is a polynomial of degree up to 
four. Since the two cases are symmetric, we discuss Di{x) in detail only. Rewrite Di{x) as 

Di{x) = diiX^ + d^x^ + d 2 x‘^ + dix + d^, 


where 


do = P-1,0 - 4 p-i,ip-i,-i, 

di = 2p_i,o(po,o - 1 ) - 4(p-i,iPo,-i +Po,iP-i,-i), 

d 2 = (po,o - 1)^ + 2pi,oP-i,o - 4(pi,ip-i,-i +pi,_ip_i,i +po,iPo,-i), 

ds = 2pi,o(po,o - 1) - 4(pi,iPo,-i +Po,iPi,-i), 

di =Pi,o -4pi,iPi,-i- 

It can be easily checked that di < 0 and da < 0. 

When Di is a polynomial of degree 4 (or d 4 / 0), there are four branch points, denoted by Xi {yi), 
i = 1,2,3,4. Without loss of generality, we assume that |a:i| < \x 2 \ < l^sl < \x/^\. When the degree of 
Di (x) is d < 4, for convenience, we let x^+k = 00 for integer k > 0 such that d + A; < 4. For example, if 
d = 3 , then X 4 = 00 . This can be justified by the following: consider the polynomial Di(x) = Di{x)/x^ 
in X, where x = 1/x. Then, 5: = 0 is a d-tuple zero of Di{x), and therefore x = 00 can be viewed as a 
d-tuple zero of Di{x). 

The following lemma characterizes the branch points of Y±[x) for all non-singular random walks, 
including the heave-tailed case, or the case of M = 0 . 

Lemma 3.1 1. For a non-singular random walk with My 7 ^ 0, Y±{x) has two branch points xi and X 2 
inside the unit cirele and another two braneh points X3 and X4 outside the unit cirele. All these branch 
points lie on the real line. More speeifically, 

( 1 ) ifPi,o > 2 y/pi,ipi,_i, then 1 < X3 < X4 < 00 ; 

( 2 ) if Pi ft = 2 ^pi,ipi,_i, then 1 < X 3 < X 4 = 00 ; 

(3) */Pi,o < 2^pi,iPi,-i, then 1 < X 3 < —X4 < 00 , where the equality holds if and only if di = da = 0. 

Similarly, 

(4) i/p- 1,0 > 2 ^p_i,ip_i,_i, then 0 < xi < X 2 < 1 ; 

( 5 ) if P-ifi = 2 ^p_i^ip_i^_i, then xi = 0 and 0 < X 2 < 1 ; 

( 6 ) if P-ifl < 2 y'p_i,ip_i,_i, then 0 < —xi < X 2 < 1 , where the equality holds if and only if di = ds = 

0 . ’ 

2. For a non-singular random walk with My = 0 (in this case M^ 7 ^ 0 since we are only considering 
the genus 1 case in this paper), either X 2 = 1 if M^ < 0; or X 3 = 1 if Mx > 0. In the latter case, the 
system is unstable. 
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Proof. We only need to prove 3. and 6. since all other proofs can be found in Fayolle, lasnogorodski 
and Malyshev [11] (Lemma 2.3.8 and Lemma 2.3.9). We provide details for 3. since 6. can be proved 
similarly. Suppose otherwise > — X4. From di < 0 and ds < 0, we obtain Di{—xz) = —d^x\ — dix^ > 
0. On the other hand, Hi(— 00 ) = —00 since ^4 < 0, which implies that Di[x) = 0 has a fifth root in 
(—oOjXs), but this is impossible. The contradiction shows that X 3 < —X 4 . It is clear that the equality 
holds if and only if di = ds = 0. □ 


Remark 3.1 Similar results hold for the branch points yi, i = 1,2,3,4, of X±{y). 

Definition 3.3 pij called X-shaped if Pij = 0 (pf^j = Oj for all i and j such that \i + j\ = 1. 

(k) 

A random walk is called X-shaped if pi ,■ and also p) f for fe = 1, 2 are all X-shaped. 

U 2,J 

Based on Lemma l3.11 we can prove the following result. 


Corollary 3.1 X 3 = —X 4 if and only if Pij is X-shaped. 


Throughout the rest of the paper, we define [x^,X 4 \ = [—oo,X4] U [x3,oo] when X4 < —1. Similarly, 
[ 2 / 311 / 4 ] = [“00,7/4] U [ 7 / 3 , 00 ]] when 7/4 < —1. We define the following cut planes: 

C-a; - C3;\[X3,X4], 

Cy = Cy \ [ 7 / 3 , ^ 4 ], 

Ca: = Ca; \ [X 3 ,X 4 ] U [xi,X 2 ], 

Cy = Cj/ \ [ 7 / 3 , 7 / 4 ] U [yi, 7 / 2 ]. 


where Ca; and Cy are the complex planes for x and 7 /, respectively. 

We now define two complex functions on the cut plane Cx based on Y±(x): 


and 


r Y.(x), if |y_(x)| < |y+(x)i, 
\ T+(x), if |y_(x)| > |y+(x)|; 

f Y+(x), if |y_(x)| < |y+(x)i, 
\ T_(x), if |y_(x)| > |y+(x)|. 


(3.2) 


(3.3) 


Obviously, Yq is the function of IT and 1+ with the smaller modulus and 1+ is the function with the 
larger modulus. 

Functions XQ{y) and Xi{y) are defined on the cut plane Cy in the same manner. 


Remark 3.2 A branch point ofY±{x) (X±{y)) is also referred to as a branch point ofYQ{x) and Yi{x) 
(Xo{y) andXi{y)). 


Remark 3.3 It is not always the case that Yq is a continuation of Y- and Yi a continuation ofYj^. 

However, for x € Ca; with a{x) / 0, To(3:) andYi{x) are still the two zeros of the kernel function h{x,y). 
Parallel comments can be made on Xq and Xi. 


A list of basic properties of Yq and li (Xq and Xi) is provided in the following lemma. 
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Lemma 3.2 1. For Ixj = 1, |Lo(x)| < 1 and |li(x)| > 1, with equality only possibly for x = ±1. For 
X = 1, we have 



for X = —1, the equality holds only if Pi,j is X-shaped, for which we have 




Yoi-l) 



2. The functions Yi{x), i = 0,1, are meromorphic in the cut plane Cx- In addition, 


(a) Yq{x) has two zeros and no poles. Hence Yq{x) is analytic in Cx', 


(b) yi(a:) has two poles and no zeros. 

(c) |lo(®)| < 1^1 (^)l; whole cut complex plane and equality takes place only on the cuts. 

3. The function 10 ( 2 ^) can become infinite at a point x if and only if, 

(a) pii = piQ = 0, in this case, x = X 4 = 00 ; or 

(b) p-ii = p -10 = 0, in this case, x = xi = 0. 

Parallel conclusions can be made for functions XQ{y) and Xi{y). 

Proof. This lemma contains results in Lemma 2.3.4 and Theorem 5.3.3 in m- First, according to 
(ii) of Theorem 5.3.3 in [H], the functions Yq and Yi defined in this paper coincide the functions Yq 
and Yi in m due to the uniqueness of the continuity. Then, all results in 1 come from Lemma 2.3.4 
and Lemma 5.3.1 in m except for the expressions for yo(~l) Yi(— 1), which can be obtained in 
the same fashion as for Yo(l) ^md Yl(1); results in 2 are given in (ii) of Theorem 5.3.3 in |11] : and the 
conclusion in 3 is the same as in (iii) Theorem 5.3.3 in m- □ 

Remark 3.4 All the above properties can be directly obtained through elementary analysis of the square 
root function. 

Throughout the rest of the paper, unless otherwise specified, we make the following assumption: 
Assumption 3 All branch points Xi and yt, i = 1,2,3,4, are distinct. 

A random walk satisfying Assumption [3] is called a genus 1 random walk. 
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Remark 3.5 This assumption is equivalent to the assumption that the Riemann surface defined by the 
kernel equation has genus 1. The Riemann surface for the random walk is either genus 1 or genus 0. 
A necessary and sufficient condition for the random walk in the quarter plane to he genus 1 is given in 
Lemma 2.3.10 in m Most of queueing application models are the case of genus 1. The genus 0 case 
can he analyzed similarly except for the heavy-tailed case, the case where M = 0. In general, analysis 
of the genus 0 case (except for the case of M = 0) could be less challenging since expressions for the 
unknown generating functions and vr 2 (y) are either explicit or less complex than for the genus 1 

case, which can immediately lead to an analytic continuation of these unknown generating functions. 
Chapter 6 of m is devoted to the genus 0 case. 

Corollary 3.2 For a non-singular genus 1 random walk, if Pij is X-shaped, then allpi^i, pi-i, P-i,i 
and are positive. 

Proof. If only one of pi,i, Pi,-i, P-i,i and is zero, then the random walk is non-singular having 

genus 0 (Lemma 2.3.10 in HI]) and if at least two of them are zero, then the random walk is singular 
(Lemma 2.3.2 in [H]). □ 

Corollary 3.3 For a stable random walk with M 0, 

1 . If p\^j is X-shaped, then p^^l and p^^^ ^ cannot he both zero; and 

(2) (2) (2) 

2. IfPij is X-shaped, then p\ [ and p_[ cannot he both zero. 

(k) 

Proof. Otherwise, Pq q = 1, /c = 1 or 2, with which the random walk cannot be stable. □ 

For our purpose, more results about functions Yq and li (Xq and Xi) are needed. Once again, we 
consider Yq and Yi. Xq and Xi can be considered in the same way. Recall that Yi, (A; = 1, 2) are defined 
on the cut plane Cx\ [x^,xfi\ U [xi,X 2 ], where two slits [xi,X 2 ] and [xs, xfi\ are removed from the complex 
plane such that the functions Y^ can stay always in one branch. Take the slit [xi,X 2 \ as an example. 
For x' G [xi,X 2 \, the limit of Yk{x) when x approaches to x' from above the real axis is different from 
the limit as x approaches to x' from below the real axis. Let x" G [xi,X 2 ] be another point satisfying 
x" > x'. By Yq\ :j^'x^' \, we denote the image contour, which is the limit of Yq{x) from above the real axis 
when X traverses from x' to x" and fro m bel ow the real axis when x conti nues to traverse back from x" 
to x'. For convenience, we say that Yq\ :^'x ' ] is the image of the contour :^'x ' , traversed from x' to x” 
along the upper edge of the slit [x', x"] and then back to x' along the lower edge of the slit. In this way, 
we can dehne the following image contours: 

■Sf = Yo[ xiX 2 ], ACext = Yolx’jxl]] (3.4) 

= ^o[ yij/2 ], -^ext = ^o[ ^3^ ], (3.5) 

respectively. Furthermore, for an arbitrary simple closed curve ^, by we denote the interior domain 
bounded by ^ and by the exterior domain. 

The properties of the above image contours provided in the following lemma are important for the 
interlace between the two unknown functions 7r(x) and 'iT 2 {y) discussed in the next section. To state the 
lemma, dehne the following determinant: 


Pll 

PlO 

Pl,-1 

Poi 

Poo 

P0,-1 

P-1,1 

P-1,0 

P-1,-1 
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Lemma 3.3 For non-singular genus 1 random walk without branch points on the unit circle, we have 
the following properties: 

1. The curve ^ and ^ext o.re simple, closed and symmetrical about the real axis in plane. 
Moreover, 

(a) If A > 0, then 

[xi,X 2 ] C C and [x3,X4] C 

(b) If A < 0, then 

[xi,X 2 ] C C and [x 3 ,X 4 ] C G^^; 

(c) If A = 0, then 

[xi,X 2 ] C = G^ and [x3,X4] C G^^. 

Entirely symmetric results hold for and .^ext- 

2. The branches Xi and Yi have the following properties: 

Yo{x) 

(a) Both XQ{y) and Yq{x) are conformal mappings: G^ — [x\,X 2 ] ^ G_^ — \yi,y 2 ]; 

Xo{y) 

(b) Xo{y) G G^ U G^^^, and Xi{y) G G'^ U G^ . Symmetrically, Yo{x) G G^ U and 

Fi(x)gG^UG^^^^; 

(c) If G^ C then 

Xq o YQ{t) = t, if t £ Gj(, 

Xo o yo(t) / t, if t£G<^j( and Xq o yo(G^) = G.^. 

Symmetrically, if Gjr C then 

YQoXQ{t) = tifteG^, 

Fo o Xo{t) / t if tGG’^ and Yq o Xo(G;^) = G^. 

Proof. A proof of the lemma can be found in Theorem 5.3.3 (i) and Corollary 5.3.5 in m- Parallel 
results when 1 is a branch point (or both 1 and -1 are branch points) can be found in Lemma 2.3.6, 
Lemma 2.3.9 and Lemma 2.3.10 of [H]. □ 

Remark 3.6 Results in this lemma can also be directly proved through elementary analysis without 
using advanced mathematical concepts used in m- 


4 Asymptotic Analysis of the Two Unknown Functions fi { x ) And 

T^2{y) 

The key idea of the kernel method is to consider all (x, y) £ B such that the right hand side of the 
fundamental form is also zero, which provides a relationship between the two unknown functions '7ri(x) 
and TT 2 {y). Then, the interlace between the unknown functions 7ri(x) and 'K 2 {y) plays the key role in the 
asymptotic analysis of these two functions, from which exact tail asymptotics of the stationary distri¬ 
bution can be determined according to asymptotic analysis of the unknown function at its singularities 
and the Tauberien-like theorem. 
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4.1 Tauberian-like theorems 


Various approaches, say probabilistic or non-probabilistic, including analytic or algebraic, are available 
for exact geometric decay. However, asymptotic analysis seems unavoidable for exact non-geometric 
decay. A Tauberian, or Tauberian-like, theorem provides a tool of connecting the asymptotic property at 
dominant singularities of an analytic function at zero and the tail property of the sequence of coefficients 
in the Taylor series of the function. In our case, an unknown generating function of a probability 
sequence is analytic at zero. Since these probabilities are unknown, in general, it cannot be verified that 
the probability sequence is (eventual) monotone, which is a required condition for applying a standard 
Tauberian theorem. The tool used in this paper is a Tauberian-like theorem, which does not require 
this monitonicity. Instead, it imposes some extra condition on analyticity of the unknown generating 
function. 

Let A{z) be analytic in Iz] < R, where R is the radius of convergence of the function A(z). We first 
consider a special case in which R is the only singularity on the circle of convergence. 


Remark 4.1 It should be noticed that for an analytic function at 0, if the coefficients of the Taylor 
expansion are all non-negative, then the radius R > 0 of convergence is a singularity of the function 
according to the well-known Pringsheim’s Theorem. 

Definition 4.1 (Definition VI.1 in Flajolet and Sedgewfick |13j l For given numbers e > 0 and 
(j) with 0 < <f < 7^12, the open domain A{(j),e) is defined by 

A{(j),e) = {z ^ C \ \z\ < 1 + e, z l,\z — 1\ > (f)} . (4.1) 

A domain is a A-domain at 1 if it is a A{(p, e) for some e > 0 and 0 < (j) < 7r/2. For a complex number 
( 0, a A-domain at ( is defined as the image C-A{(p,e) of a A-domain A{4>,e) at 1 under the mapping 

z C,z. A function is called A-analytic if it is analytic in some A-domain. 

Remark 4.2 The region A{4>,£) is an intended disk with the radius of 1 -\- e. Readers may refer to 
Figure VI. 6 in IWf for a picture of the region. Throughout the paper, without otherwise stated, the limit 
of a A-analytic function is always taken in the A-domain. 


Theorem 4.1 (Tauberian-like theorem for single singularity) Let A{z) = be ana¬ 

lytic at 0 with R the radius of convergence. Suppose that R is a singularity of A{z) on the circle of con¬ 
vergence such that A{z) can be continued to a A-domain at R. If for a real number a ^ {0, —1, —2,...}, 


lim (1 

z^R 


zlRTA{z) = g, 


where g is a non-zero constant, then. 


r(a) 

where r(Q:) is the value of the gamma function at a. 




Proof. This is a immediate consequence of Corollary VI. 1 in [13] after the transform z Rz. □ 

For the random walks studied in this paper, we will prove that the unknown generating function 
T^i{x) {'K 2 {y)) has only one singularity on the circle of its convergence, except the X-shaped random 
walk for which the convergent radius R and —R are the only singularities. To deal with the later case, 
we introduce the following Tauberian-like theorem for the case of multiple singularities. 
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Theorem 4.2 (Tauberian-like theorem for multiple singularities) Let A{z) = be an¬ 

alytic when \z\ < R and have a finite number of singularities fk, h = 1,2,... ,m on the circle \z\ = R of 
convergence. Assume that there exists a A-domain Aq at 1 such that A can be continued to intersection 
of the A-domains Cfc at fk, k = 1,2 ,... , m: 


^ = nr=i(a-Ao). 


If for each k, there exists a real number Ok ^ {0, —1, —2,...} such that 


lim (1 - z/CkT^A{z) = 9k, 


where gu is a non-zero constant, then, 


E 


9k 


^ r(afc) 


n 




Proof. This is an immediate corollary of Theorem VI. 5 in [T3] for the case where Ok is real, fik = 0, 
c^kiz) = Tk{z) = (1 - 2 ;)““'' and ak,n = □ 

4.2 Interlace of the two unknown functions 'Ki{x) and TV 2 {y) 

The interlace of the unknown functions t^i{x) and TT 2 {y) is a key for asymptotic analysis of these 
functions. Let 

Pq = {x e C :|x| = a}, 

Da = {x : \x\ < o}, 

Da = {x : |x| < a}. 

When a = 1, we write P = Pi, D = Di and D = Di. 

We fist state two literature results on the continuation of the functions tti{x) and 'K 2 {y). 

Lemma 4.1 (Theorem 3.2.3 in |11| I For a stable non-singular random walk having genus 1, 'Ki{x) 
is a meromorphic function in the complex cut plane Cx- Similarly, TT 2 {y) is a meromorphic function in 
the complex cut plane Cy. 

This continuation result is crucial for tail asymptotic analysis. The following intuition might be 
helpful to see why such a continuation exist. When the right hand side of the fundamental form is zero, 
the X and y are related, say through the function Yq{x). Therefore, x^ is the dominant singularity if 
there are no other singularities exist between (1,X3). Based on the expression for irfix) obtained from 
the fundamental form, all other singularities come from the zeros of hi{x,YQ{x)), which are poles of 
7ri(x), or the singularities of Tr 2 {YQ(x)). A similar intuition holds for the function T^ 2 {y)- Based on the 
above intuition, it is reasonable to expect Lemma l4.ll 

Remark 4.3 An analytic continuation can be achieved through various methods. In 111]/ and |j5]/ . it 
was proved in terms of properties of Riemann surfaces. In JWf and IWj . direct methods were used for a 
convergent region. For some cases, a simple proof exists by using the property of the conformal mapping 
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Yq or Xq. For example, for the case of My > 0 and Mx < 0, we know, from Lemma W^ l. Lemma WfT\ 
and Lemma WM -2 respectively, that |5^o(2^)| < 1 for jxl = 1, > 1 and Yq is analytic in the cut plan. 

Therefore, it is not difficult to see that we can find an e > 0 such that for |x| < 1 + e, the function 
'^ 2 {Yo{x)) in is analytic, which leads to the continuation ofni^x). 


Lemma 4.2 (Lemma 2.2.1 in [llj i Assume that the random walk is ergodic with M 7 ^ 0 and the 
polynomial h{x,y) is irreducible. Then, exists an e > 0 such that the functions vri(x) and 7r2(y) can be 
analytically continued up to the circle Fi+g in their respective complex plane. Moreover, they satisfy the 
following equation in -Df+g H B: 

hi{x, y)7ri(x) + h 2 {x, y)TT 2 {y) + hQ{x, y)7ro,o = 0. 


Proof. The analytic continuation is a direct consequence of Lemma l4.ll and the equation is directly 
from the fundamental form. □ 


Theorem 4.3 1. Function Tr 2 {YQ{x)) is meromorphic in the cut complex plane C^. Moreover, ifYo^x^) 
is not a pole of 'K 2 {y), then xq is x^om of tt 2 {Yq{x)) and there exist e > 0 and 0 < f < 'njl such that 

lim 7r2(yo(a:^)) = 'n' 2 {Yo{x 3 )) and lim -k' 2 {Yq{x)) = 7r2(To(a;3))- 


Similarly, 7ri(Xo(y)) is meromorphic in the cut complex plane Cy. Moreover, i/Xo(y 3 ) is not a pole 
of TTi{x), then ys is ydom of TTi^Xo^y)) and there exist e > 0 and 0 < < 7 r /2 such that 

lim 7 ri(Xo(y)) = 7ri(Xo(y3)) and lim 7 r[{Xo{y)) = TT[{Xo{y 3 )). 
y^ys y^ys 


2. In cut plane Cx, equation 

hi{x,Yo{x))TTi{x) + h2{x,Yo{x))7r2{Yo{x)) + ho{x,Yo{x))7rofi = 0 

holds except at a pole (if there is any) of 7ri{x) or 7r2{YQ{x)). Therefore, 

-h2{x,Yo{x))7r2{Yo{x)) - ho{x,YQ{x))7:ofi 
hi{x,Yo{x)) 

except at zero of hi{x,YQ{x)), or at a pole (if there is any) of iti{x) or it 2 {Yq{x)). 
Similarly, in the cut plane Cy, equation 

hi{Xo{y),y)Tri{Xo{y) + h2{Xo{y),y)Tr2{y) + ho(Xo(y), y)7ro,o = 0 

holds except at a pole (if there is any) of 7r2{y) or 7ri(Xo(y)). Therefore, 

-hi{Xo{y),y)7i:i{Xo{y)) - ho{XQ{y),y)TTofl 

=-mwT-• 

except at a zero 0 //i 2 (Xo(y), y), or at a pole (if there is any) of n 2 {y) or 7ri(Xo(y)). 


(4.2) 


(4.3) 


(4.4) 


(4.5) 
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Proof. We only prove the result for functions of x and the result for functions of y can be proved in 
the same fashion. 

1 . From Lemma ED and Lemma HU Yq{x) is analytic in the cut complex plane C^; and 7r2(y) is 

meromorphic in the cut complex plane Cy, which implies T^ 2 iXo{.x)) is meromorphic in if Lo(3^) ^ 
[?/3,?/4]. According to Lemma [3^ 2fbi. for all x S C^;, Yq{x) G and according to Lemma 

1) [ys^Vi] C {G_^ U which conhrms L()(2^) ^ [ys^yi]- From the above, we have 'K 2 {y) is analytic 

at 10(3^3)) then the limits in 1 . are immediate results of the analytic properties of '7r2(lo(3^))- 

2. Since both 7ri(x) and 7r2(Lo(3;)) are meromorphic (proved in 1.) and Yq[x) is analytic fLemma l.I.ip 

in Cx, equation (14.21) in the cut plane Ca, except at the poles of 7ri(x) or 7r2(Lo(3^))- D 

Remark 4.4 Let us extend the definition of Tri{x) to x = x^ by vri(x3) = lim3;^a;3 7ri(a:) for x in the 
cut plane. We say that X3 is a pole if the limit of 7ri{x) is infinite as x ^ X3 in the cut plane. 

According to the above interlacing property and the Tauberian-like theorem, for exact tail asymp¬ 
totics of the boundary probabilities vr^^o and we only need to carry out an asymptotic analysis at 
the dominant singularities of the functions vri(x) and vr 2 (y), respectively. There are only two possible 
types of singularities, poles or branch points. We need to answer the following questions: 

Ql. How many singularities on the circle of convergence (dominant singularities)? 

Q2. What is the multiplicity of a pole? 

Q3. Is the branch point also a pole? 

For the random walk considered in this paper, we will answer all these questions. We will see that 
on the convergent circle, there is only one singularity or there are exactly two singularities. For the 
former, Theorem 14.11 will be applied, and for the latter. Theorem 14.21 will be applied. 

4.3 Poles of 7ri(x) 

Farallel properties about poles of the function TT2{y) can be obtained in the same fashion, which will 
not be detailed here. 

Lemma 4.3 1. Let x & G^r\{Dy, then the possible poles of 'Ki{x) inG^r\{Dy are necessarily zeros 
of hi{x,Yo{x)), and |To(a;)| < 1 - 

2 . Let y G Gjr r\ {Dfi, then the possible poles of iT2{y) in G^ n {Dfi are necessarily zeros of 
h2{Xo{x),y), and |Xo(y)| < 1. 

Froof. 1. When x G then Yq{x) = y G [yi,y 2 ]- From Lemma EH for |x| = 1, |Fo(a;)| < 1- For 
X G G_^r\ {Dfi, it follows from the maximum modulus principle, we have |Lo(®)| < 1- Hence, tt 2 {Yq{x)) 
is analytic in G^^ n {Dfi. From Theorem 14.31 if hi{x,YQ{x)) 0, equation (14.31) holds, which implies 
that the possible poles of 7ri(x) in G^ (H [Dfi are necessarily zeros of hi{x,YQ{x)). 

2. The proof is similar. □ 

Theorem 4.4 Let Xp be a pole of tti(x) with the smallest modulus. Assume that \xp\ < X 3 . Then, one 
of the follow two cases must hold: 
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1. Xp is a zero of hi {x,Yq{x)); 

2. yo = Yo{xp) is a zero of h 2 {Xo{y),y) and \yo\ > 1. 

Parallel results hold for a pole of Tr 2 {y). 

Proof. Suppose that Xp is not a zero of hi{x,Yo{x)). According to equation (j4..S[) in Theorem 14.31 Xp 
must be a pole of 'K 2 {Yq{x)) and |yo| > 1- Furthermore, by Lemma 1431 Xp ^ G^. If yo is not a zero 
of h 2 (Xo(y),y), according to equation (14.5j) in Theorem 14.31 yo must be a pole of 7ri(Xo(y)), that is, 
xq = Xniftfo) is a pole of 7ri(x). It follows from Lemma 031 that xq = Ao(yo) is a zero of hi{x,YQ{x)) 
if Xq € Gj^. There are two possible cases: A > 0 or A < 0. If A > 0, by Lemma [TSl lfai and 2(c), 
Xq e G^. In the case of A < 0, according to Lemma 1531 llbL 1(c) and 2(b), we also have xq G G^. 
However, this case is not possible, since otherwise according to Lemma 13.31 1 we would have xq = Xp or 
Xq = —Xp, both leading to a contradiction. This completes the proof. □ 

Remark 4.5 IFe will show in the next subsection that a pole of tti{x) with the smallest modulus in the 
disk |x| < X3 is real. 

4.4 Zeros of/ii(a;, 10 ( 2 ;)) 

In this subsection, we provide properties on the zeros of the function hi{x,YQ{x)). The main result is 
stated in the following theorem. 

Theorem 4.5 For a non-singular random walk having genus 1, consider the following two possible 
cases: 

1. Either pi^j or p[^j is not X-shaped. In this case, either hi{x,YQ{x)) has no zeros with modulus in 
( 1 ,X 3 ], or it has only one simple zero, say x*, with modulus in ( 1 ,X 3 ], and x* is positive. 

2. Both Pij and p\^j are X-shaped. In this case, either hi{x,Yo{x)) has no zeros with modulus in 
( 1 , 0 : 3 ], or it has exact two simple zeros, namely, 0 :* > 0 (with modulus in ( 1 , 0 : 3 ]^ and —x*, both are 
zeros of hi{x,Y q{x)) or both are zeros of a{x)hi{x,Yi{x)). 

With this theorem and Theorem 14.41 we are able to apply the Tauberian-like theorem to characterize 
the tail asymptotic properties for the boundary probability sequence TTnfl- To show the above Theorem, 
we need the following several lemmas and two propositions. Instead of directly considering the function 
/o(x) = hi{x,Y(i{x)), we consider a polynomial f(x}, which is essentially the product of fo(x) and 
fi(x) = hi(x,Yi(x)): 

f(x) = Mx}fi(x}, 

where fi(x) = a(x)fi(x). It is easy to verify, by noticing 

Yo(x)Yi(x) = and yo(a:) + Ti(x) = 

a{x) a{x) 

that 

f(x) = a(x)bf(x) — b(x)bi(x)ai(x) + c(x)af(x) (4.6) 

= dex^ + d^x^ Y d4x‘^ + dsx^ + d2X^ Y dix-h dg. (4.7) 
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Hence, a zero of fi{x), z = 0,1, has to be a zero of /(x), and any zero of f{x) is either a zero of fo{x) 
or a zero of fi{x) = a{x)fi{x). 


We can also write 

where 


with 


/(x) = a(x)[oi(x)]^R_(x)it’+(x), 

(4.8) 

RU^) - Fix) ± 

(4.9) 

MX) Hx) 

oi(xj 2a(xj 

(4.10) 


Remark 4.6 1. It can be easily seen that both fo{x) and fi{x) are analytic on the cut complex plan. 
In fact, the analyticity of fo{x) is obvious and the analyticity of fi{x) is due to the cancellation of the 
zeros of a{x) and the pole of fi{x). 


All proofs for Lemmas I4.4H4.7I and for Proposition 14.11 and Proposition 14.21 are organized into Ap¬ 
pendix lAl 


Lemma 4.4 1. (a) yo'(l) = ^ if My < 0; (b) Y({1) = ^ zf My > 0; and (c) ^(1) = ^ 0 ( 1 ) 
and X = 1 is a branch point ofYi{x) and Yo{x) if My = 0. In this case, Y{{1) and ^ 0 ( 1 ) do not exist. 
Parallel results hold for functions Xk{y). 

2. If My / 0, then f{x) has at least one non-unit zero in [x 2 ,xf\ and 1 is a simple zero of f{x). 
Parallel results holds for the case of M^ 7 ^ 0. 


Lemma 4.5 1. Let z be a branch point ofY^fx). If f{z) = 0, then z cannot be a repeated root of 
f{x) = 0 . 

2. f{x) (therefore both fo{x) and fi{x)) has (have) no zeros on the cuts, except possibly at a branch 
point. More specifically, f{x) < 0 if a{x) < 0 and f{x) > 0 if a{x) > 0. 

3. fo{x) and fi{x) have no common zeros except possibly at a branch point or at zero. 


4 . Consider the random walk in Theorem 4.5-1. If fo{x) has a zero in [—x^,- 
additional (different) zero in [—X 3 ,— 1 ). 


-1), then /o(x) has an 


5. For the random walk in Theorem \4.5[ l, if \x\ G (IjXs], then |lo(~|3;|)| < yo(|3;|)- 


Lemma 4.6 Consider the random walk in Theorem \4.5\ -l. If My < 0, then x = 1 is the only zero of 
fo{x) = hi{x,YQ{x)) on the unit circle |x| = 1. If My > 0, then fo{x) has no zero on unit circle |x| = 1. 


Remark 4.7 From the proof of Lemma 4-6. we can see that for the random walk considered in Theo- 
rem\4.5[2, fo{x) has no zeros with non-zero imaginary part on the unit circle. 


The proof of Theorem 14.51 is based on detailed properties of the function /(x) and also the powerful 
continuity argument to connect an arbitrary random walk to a simpler one. For using this continuity 
argument, we consider the following special random walk. 
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Special Random Walk. This is the random walk for which pij is cross-shaped (or pij = 0 
whenever \ij\ = 1), and ^ o “ prove the counterpart result to Theorem 14.51 for the 

Special Random Walk. 


Proposition 4.1 For the Special Random Walk, the following results hold: 

1. f{x) = 0 has six real roots with exact one non-unit root in [ 3 : 2 , 3 : 3 ]. More specifically, two roots 
are zero, two in {x 2 X^\, one in (—oo,xi], and one in [x 4 ,oo). 

2. If fo{x) has a zero, say x*, in (1, xfi\, then x* is the only zero of fo{x) with modulus in (1, X 3 ]. 
Furthermore, fo{x) has no other zeros with modulus greater than 1 except possibly at x = x^. 


For the random walk considered in Theorem 14.51 2. we hrst prove the following results. 

.( 1 ) 


Lemma 4.7 For the random walk considered in Theorem 4-5-2 (or both pij and p\ ■ are X-shaped), 
/(I) = /(—I) = 0 , and f{x) = 0 has two more real roots, say 0 < xq 7 ^ 1 and —xq, and two complex 
roots. 


Proposition 4.2 For the random walk considered in Theorem \4.5[ 2 (or bothpij andp^^j are X-shaped), 
either the two complex zeros of f{x) are zeros of fi{x) = a{x)fi{x) or they are inside the unit circle. 


Proof, of Theorem 14.51 1. For the random walk considered here (either pij or p\^j is not X-shaped), 

P = {p-l-l,P 0 -l, Pi,- 1 , P- 1 , 0 , P 0 fi,P 0 ,l,P-l,Pi, 1 ), 

„(1) _ /XI) XI) XI) XI) XI) Xi)^ 

P ~ W-i,o^Po,o^Po,i^P-i,i^Po,i^Pi,i)- 


let 


Define 


A = {(p,P^^^) : 0 < Pij,p^ij < 1 and ^Pi,j = J^PiJ = 


*>3 


U3 


For an arbitrary random walk for which either pij or p^j is not X-shaped, let p be the corresponding 
point in A. We assume that My < 0 for the random walk p (and a similar proof can be found for the 
case of My > 0). Let po be an arbitrarily chosen point in A corresponding the Special Random Walk. 
We prove the result by contradiction. Suppose otherwise that the statement were not true. There would 
be three possible cases: (i) Im(x*) 0; (ii) —X 3 < x* < —1; and (hi) there exists xq G ( 1 , 3 ^ 3 ] with 
xq / X* such that /o(xo) = 0. 

Case (i). Clearly, x* is also a root of f{x) = 0. Choose a simple connected path ^ in X to connect 
p to po such that on I (excluding p, but including po) My < 0. The zeros of f{x) as a function of 
parameters in A are continues on I. There are two possible cases: (a) the zero function xq{9) (with 
xo(p) = X*) never passes the unit circle when 6 travels from p to po; and (b) xo{6) passes the unit circle 
at some point 6 ^ i. 

If (a) occurs, let 9o be the^first point at which xo{9) = xo{9), where xo{9) is the zero function with 
xo(p) = X*. If X* is a zero of /i, then /o and /i would have a common zero xo(9o) = x(9o) at 9o, which 
contradicts Lemma 14.51 3. Hence, the only possibility is that x* is also a zero of /q. From 9o on, both 
xo(9) and x(9) shouldjilways be zeros of fo, since otherwise only at a branch point a zero of fo could be 
switched to a zero of /i and all branch points are real, which means that xo(9) = x(9) is a branch point 
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and a multiple roots, contradicting to Lemma 03]- 1. As Oq approaches po, it leads to a contradiction 
that two zeros of /o are in (IjXsj. 

If (b) occurs, we can assume that when xo{6) pases the unit circle it is a zero of /o based on the 
proof in (a). Then, /o has two zeros since 1 is always a zero of /o independent of the parameters (or 9) 
when My < 0, which is a different zero from xq{9). This contradicts to the fact that /o has only one 
zero at the unit circle. 

Case (ii). In this case, fo{x) would have another zero in [— X 3 , —1) at p according to Lemma [45l 4. 
Consider the same two cases (a) and (b) as in (i). We can then follow a similar proof to show that case 
(ii) is impossible. 

Case (iii). A similar proof will show that the case is impossible. 

2. This is a direct consequence of Lemma 14.71 and Proposition 14.21 □ 

The following Lemma gives a necessary and sufficient condition under which /o(x) = hi{x,YQ{x)) 
has a zero in ( 1 , X 3 ]. 

Lemma 4.8 Assume My 7 ^ 0. We have following results: 

1. ///o(x 3 ) > 0, /o(x) has a zero in ( 1 ,X 3 ]; 

If foixs) < 0, fo{x) has no zeros in ( 1 ,X 3 ]. 

Proof. 1. There are two cases: My > 0 or My < 0. If My > 0, then /o(l) < 0, which leads to the 
conclusion. If My < 0, then /o(l) < 0, which also leads to the conclusion since /o(l) = 0 and fo{x^) > 0. 

2. Again there are two cases; My > 0 or My < 0. By simple calculus, in either case, we obtain 

that if fo{x) = 0 had a root in ( 1 ,X 3 ], then it would have another root in ( 1 , 0 : 3 ] since foixs) < 0. This 
contradicts to Theorem 14.51 □ 

4.5 Zeros of h 2 {Xo{y),y) 

Following the same argument in the previous subsection, we have the following result: 

Theorem 4.6 For a non-singular random walk having genus 1, consider the following two possible 
cases: 

i‘2') 

1. Either pi^j or pi ■ is not X-shaped. In this case, either /i 2 (Ao(y),y) has no zeros with modulus 
in ( 1 , 2 / 3 ], or it has only one simple zero, say y*, with modulus in ( 1 , 2 / 3 ], and y* is positive. 

( 2 ) 

2 . Both pij and pl ■ are X-shaped. In this case, either h 2 {XQ{y),y) has no zeros with modulus in 
( 1 , 2 / 3 ], or it has exact two simple zeros, namely, y* > Q (with modulus in {l,y^]) and —y*, both are 
zeros of go{y) or both are zeros of gi{y), where 

90 {y) = h 2 {Xo ( 2 /), 2 /) and gi (y) = /12 {Xi {y),y). 

From the above analysis, we know that if hi{x, Yo(x)) has a zero in (1, xf\, then such a zero is unique. 
Similarly, if /i 2 (Ao( 2 /), 2 /) has a zero in ( 1 , 2 / 3 ], then such a zero is unique. For convenience, we make the 
following convention: 

Convention 1 Let x* be the unique zero in ( 1 , 0 : 3 ] of the function hi{x,YQ{x)), if such a zero exists, 
otherwise let x* = 00 . Similarly Let y* be the unique zero in (1, 2 / 3 ] of the function h 2 {XQ{y),y) if such 
a zero exists, otherwise let y* = 00 . 
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According to Theorem 14.41 the unique pole in ( 1 , 3 : 3 ] of 7ri(x) is either x*, or the image of the pole 
under Yq is a zero of h2iXo{y),y). Our focus in this subsection is on this special case of y*. 

Theorem 4.7 If the pole in (1,X3] of Tri(x) is not x*, then, it, denoted by xi, satisfies: 

1 . xi =Xi{y*), where y* is the unique zero in ( 1 , 2 / 3 ] of the function h2{XQ{y),y); 

~ f2') 

2. xi is the only pole of Tri{x) with modulus in ( 1 , 2 / 3 ], exeept for the case where both pij and p\ ■ 

are X-shaped, for whieh —xi is the other pole of tti{x) with modulus in ( 1 , 2 / 3 ]- 

Proof. 1. Let x be the solution of y* = To(x). Then, x = xq = Xo{y*) 01 x = xi = Xi{y*). If 
y* then x = xq so that y* = lo(Alo(2/*))- this case, by Lemma [T3l xq < 1. If y* £ then 

5: = xi so that y* = YQ{Xi{y*) and xi G 

2. It follows from the fact that the zero, y*, of h2{XQ{y),y) in ( 1 , 2 / 3 ] is unique and the fact that 
y* = Yq{x) has only two possible solutions xq < \ and xi. In the case where both pij and pi^ are 
X-shaped, —y* is the other zero of /i 2 (Ao( 2 /), 2/) with either —y* = Yq{—xi) or —y* = Yq{—xq). □ 

Corollary 4.1 Let x be a solution of y* = Yq{x). In order for x to be in (1,X3] we need y* G 
Furthermore, we have y* < ys- 

Proof. The first conclusion is directly from the proof to Theorem 14.71 and the second one follows from 
that fact that by Lemma [33I l and Lemma [T3l 2fbL there exists no a: G (1, X 3 ] such that y* = ys = Yq{x). 
Therefore, we should have y* < ys- □ 

Convention 2 Let xi = Xi{y*) if the unique zero y* in ( 1 , 2 / 3 ] of the function h2{XQ{y),y) exists, 
otherwise let xi = 00. 


4.6 Asymptotics behaviour of 'Ki{x) and 'K 2 {y) 

In this subsection, we provide asymptotic behaviour of two unknown functions vri(a:) and TT 2 {y)- We 
only provide details for 7 ri(a:), since the behaviour for TT 2 {y) can be characterized in the same fashion. 

It follows from the discussion so far that: 

( 1 ) If pij is not X-shaped, then, independent of the properties of p\^j and , there is only one 

dominant singularity, which is the smallest one of x*, xi and X 3 . Here x*, xi and X 3 are not 
necessarily all different. 

( 2 ) If pij is X-shaped, then both X 3 and —X 3 are branch points. 

(a) If p\^j is not X-shaped, then hi{x,YQ{x) has either no zero or one zero x* in ( 1 ,^ 3 ]; and if 
p\^j is X-shaped, then hi{x,YQ{x) has either no zero or two zeros x* G ( 1 , 3 : 3 ] and —x*. 

(b) Similar to (a), h 2 {XQ{y),y) has either no zero in ( 1 , 2 / 3 ] or one zero y* in it. For the latter, if 

i‘2') ~ 

Pi ■ is not X-shaped, then xi = Xi{y*) is the only pole of 7 r 2 (Tb( 3 :)) with modulus in ( 1 ,^ 3 ]; 

and if pY is X-shaped, then xi = Xi{y*) G ( 1 , 3 : 3 ] and — 3:1 = Xi{—y*) are the only two 
poles of 7r2(To(3:)) with modulus in ( 1 , 3 : 3 ]. 
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Therefore, in case (2), we either have only one dominant singularity or exactly two dominant 
singularities depending on which of x*, xi and X3 is smallest and the property of A: = 1, 2. 


The theorem in this subsection provides detailed asymptotic properties at a dominant singularity 
for all possible cases. Let Xdom be a dominant singularity of 7ri(x). Clearly, |xrfom| = x*^ \xdom\ = xi 
or I Xdom I = X 3 . To state this theorem for the cases where Xdom = notice that through simple 

calculation we can write 

hi{x,Yo{x)) =pi{x) + qi{x)./l -(4.11) 

V ^dom 

Yq{x) = p{x) + q{x)./l - (4.12) 

V ^dom 

Yoixdom) - Yo(x) = (l - p*{x) - q{x)Jl -(4.13) 

\ ^dom J V ^dom 

hi{x,Yo{x)) - hi{xdom,Yo{xdom)) = (1 - —] pI{x) + qi{x). I -(4.14) 

V ^dom J V ^dom 

where 


p{x) 

pI{x) 


-bjx) 

2a(x) 


^dom 




p*{x) = 


2 a,(^x') 2 cl(^x 


1 


^do7 


■(x Xdom) 


. , -b(x)ai(x) , , , 
Pdx) = — \ \ , +bi{x), 


2a(x) 


fli(x) Oi(Xdom) T (x) fel(33dom)\ 

^dom ^ / 


q{x) 



/ Di{x) 

2a(x) 

1/ 1 

V ^dom 

1 

Di{x) 

2 a(x) Y 

^dom 


if Xdom a;3, 
if Xdom a^S) 


and qi{x) = ai{x)q{x). 
Define 


L(x) 

L{y) 


[/i2(x,yo(x))7r2(yo(x)) + ho(x,yo(x))7ro^o]hi(x, Yi(x))a(x) 

x/'(x) 

[hi{Xo{y),y)7i:i{Xo{y)) + ho{Xo{y),y)7rofi]h2{Xi{y),y)a{y) 

yg'iy) 


where/(x) = a(x)/ii(x, To(a;))^i(a;, yL(x)) is a polynomial defined in Section im and aiv) = a{y)h 2 {Xi{y),y) 
h 2 {XQ{y),y) is the counterpart polynomial for function /12. 

The following Theorem shows the behaviour of 7ri(x) at Xdom- Recall that yo = ^0(2^*)- 


Theorem 4.8 Assumed that both h 2 {x*,YQ{x*))7r2{YQ{x*)) + /io(x*, lo(y*))7ro,o / 0 and /ii(Xo(yo)) yo) 
7r(Xo(yo)) + b^o{Xo{yo),yo)TTofl / 0. For the function vri(x), a total of four types of asymptotics exist 
as X approaches to a dominant singularity of tti{x), based on the detailed property of the dominant 
singularity. 


Case 1: If \xdom\ = x* < min{xi,X3}, or \xdom\ = xi < min{x*,X3}, or \xdom\ = x* = xi = X3, 


then 


lim I 1 — 

^ ^^dom 


7ri(x) — ^ 


^dom 
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where 




' L{xdom), _ i/x* < min{xi,X3}; 

-h2{xdom,yo)yoL{yo) 


hi {xdom: yo')^(){xdoin')Xdom 

h2{xdomiyo)L{yo)yo 


' Qlixdom)Qixdom) 

with yo = Yo{xdom)- 

Case 2: If \xdom\ = x* = xsKxi or \xdom\ =xi= X 3 < x*, then 


, if xi <mm{x*,X 3 }-, 
if X* = XI = X3, 


lim x/xdom'^lix) — Cq 2iXdom) 1 

— - >1 


where 


C0,2(^dom) — ^ 


^2(^dom) yo)^2( 2 / 0 ) “1“ h()(^Xdomi yo)^0,0 
~fli {Xdom) 
h2{xdom,m)yoL{yQ) 

, hl{Xdomiyo)y.{Xdom) 


, ifx*=X3<Xi- 
ifxi = X3 < X*, 


with yo = Yo{xdom)- 

Case 3: If\xdom\ = X 3 < min{5:i, x*}, then 


lim i/l- x/XdomT^lix) = Co^siXdom), 

^ ^^dom 


where 


C0,3(^c/om) — 


q{Xdom)xlom ^ 


2 dy 

Case 4' If la^doml = x* = xi < X 3 , then 


lim I 1 — 

X yx^dor 


h2{Xdom,y)'^2{y) + ho{Xdom,y)'^{),0 
h\ {Xdora'i 2/) 


7 ri(x) — C{)^i\{Xdom) 1 


y—^0 {x^dom) 


Xdom 


where 


CQ, 4 (^^ 0771 ) 


h2{xdomiyo)[hi{xo,yo)T^i{xo) + ho{xo,yo)]7ro,o] 


x*'^h[{xdom, yo)Yo{xdom)h 2 {XQ{yo),yQ) 


with yo = Yo(xdom) and xq = ^o(yo)- 


Proof. Case 1. If x* < xi, then Xdom is not a pole of 7 r 2 {Yo{x)). According to Theorem 14.51 Xdom is 
a simple pole of '/ri(x). Prom equation (|4.3n in Theorem 14.31 and Lemmas 14.41 and 14.51 we have 


7ri(x) = 


-/j,i(x,lo(x)) 7 r 2 (Yb(x)) - ho{x,Yo{x))TTofl 
hi{x,Yo{x)) 

-[hi(x,yo(x)) 7 r 2 (yo(x)) + ho(x,yo(x)) 7 ro,o]hi(x, Yi)a(x) 

fix) 

-[hi{x, Yo{x))TT 2 iYoix)) + ho(x, yo(a;))7ro,o]/ii(x, Yi)a(x) 


(x - Xdom)f*ix) 


where f*ixdom) = f'ixdom) + 0. It follows that 


lim I 1 — 

X yXflor 


X 


Xdom 


7 ri(x) = L{xdom)- 
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Similarly, if xi < x*, following the same argument used in the above, we have T^ 2 {y) = 

L(yo) and 


lim I 1- 

V XdQui 


7ri(x) 


= lim 

^ ^^dom 


-h2{x,Yo{x)) (^1 - 7r2{Yo{x)) - (1 - ^^)ho{x,Yo{x))7rofl 


I ypfr) 

j^^^hi{x,Yo{x)) 

^dom 


^ -h2{xdom,yo)L{yo)yo 

^lixdom} yo)YQ(^x) {xdom)Xdom 

In the case of x* = xi = X 3 = \xdora\, we first have lima,^^;^^^ (l - 'n' 2 iYo{x)) = L{yo). Then, 

using equations ()4.13l) . ()4.14l) and the expression for hi{x 3 ,yo), we obtain 


lim I 1 - 

^^^dom \ ^dom 


7ri(x) 

/ 


= lim 

^ ^^dom 


-h 2 {x,Yo{x)) \_^ [(1 - 7 r2{Yoix))] - J. 

VQ ^ 


- Jl- 


^dom 


/io(x,yo(a^))7ro,o 


hi{x,YQ{x))/ 


^don 


L{yo)h2ixdom,yo)yo 

Ql {xdom)y{xdom) 


Case 2. If x* = xs^ then hi{xdom,yo) = Pi{xdom) = 0, using equations (I4.3p . (|4.1ip . p4.I3p and 
p4.I4p . we can rewrite vri(x) as 


vri(x) = 


-h2ix,Yo{x))7r2{Yo{x)) - ho{x,Yo{x))7rofl 


yT^^xTx^ ^I - x/xdomP*i(x) + qi (x) 


It follows that 


A-- , X -h2(x,yo(a;))vr2(yo(a;)) - ho(x,yo(a;))7ro,o , ^ 

lim VI - x/xdom 7 ri(x) = ^ lim - p - - 1 -^ = co, 2 (a:dom)- 

Vl - xlxdomP\{x) + g'i(x) 


X^Xd. 


X^Xdo 


Note that qi{xdom) V 0. 

Similarly, if xi = X3, then yo is a pole of 712 ( 1 /), which gives limy^y^ fl - n{y) = L{yo). Again, 
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using equations (14.3p . ()4.1ip . (|4.13l) and (|4.14l) . we obtain 


lim , /1 — 

X yX(io', 


X 


^dom 


- 7 ri(x) 


= lim 


-h2{x,Yo{x)) ^l^ (l - TT^iYoix)) - ^l-^hoix,Yoix))7ro,o 


hi{x,Yo{x)) 


lim 


yo^n - 


^do 


_ h2{Xdom,yo)L{yo) 

hl{Xdom,yo) (1 - x/xdom)P*{x) - q{x)^l - xjxdom 

^ h 2 {xdom,yo)L{yo)yo 
{xdom 1 yo)q{xdom) 

Case 3. Let 


T{x,y) = 


h2{x, y)'i^2{y) + hQ{x, y)'Kofl 
-hi{x,y) 


Then, 


with 


,, , 5r , dTdYoix) 
= p'{x) + q'{x)y^l - xjxdora “ 


q{x) 


and 


dT 
dy 


‘^Xdom "x/l X jXdom 

dT ^ d2{y)'^2{y) + ao{y) + [a'iix)y + b[{x)]T{x,y) 
dx -hi{x,y) 

+ h^ix, + *^r(x, y) 


-hi{x,y) 

where p{x) and q{x) are defined by equation (|4.12p . Since limj;^^:^;^^ a/1 - x/xdom^ = 0, lima;.^^;^^^ 
xjxdon f^l^"' = - ‘^txalZ^ continuous at {xdom,Yo{xdom)), 


hm ^1-x/xdomT^'lix) = -^^7^^l(-3,yo) 

^'^dom dom 

_ ^(^dom) dTiXfiQYa'i !?/) i 


—j/0 — C0,3(^dom)* 


(4.15) 

2xdom dy (4.16) 

It is easy to see c^fiixdom) / 0; since otherwise Tr'iixdom) < oo, which contradicts the fact that x^ is a 
branch point of vri(x). 

Case 4. From equation (14.31) and (|4.5I) in Theorem 14.31 we have 

h2{x,Yo{x))hi{Xo{Yo{x)),Yo{x))7ri{Xo{Yo{x))) + N{x) 


7ri(x) = 


hi{x,Yo{x))h2{Xo{Yo{x)),Yo{x)) 


where 


Since 


N{x) = [h 2 {x,Yo{x))ho{Xo(Yo{x)),Yo{x)) -/i2(^o(4o(a;)), lo(x))/io(x, lo(a;))]vro,o- 

lim fei(x.y„(x)) - ^ j 


^ ^^dom X ^dOTTL ^ ^^do 


^ ^dom 
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and 


h2{Xo{Yo{x)),Yo{x)) h2{Xo{Yo{x)),Yo{x))-h2{Xo{yo),m) 

lim - = iim - 


^ ^^dom X X(^Q'fYi 


^ ^^dom X XdQ'fji 

^ (^rfom)^2 (^0 (?/o)) ?/o)) 


we obtain 


lim 


1 - 


^do 


x^xdom hi{x,Yo{x))h2{Xo{Yo{x)),Yo{x)) xl^^h[{xdom,yo)YQ{xdom)h2{Xo{yo),yo)' 
which yields 


lim 1- 

V Xdom ) 


's 


^l(^) — CQ^iixdom)- 


□ 


Remark 4.8 It should be noted that the above theorem provides the asymptotic behaviour at a dominant 
singularity, either positive or negative. 

Corollary 4.2 If h 2 {x*,Yo{x*))'K 2 {Yo{x*))+hQ{x*,YQ{y*))-Kofl = 0 or hi{Xo{yo),yo)'K{Xo{yo))+ho{Xo{yo), 
yo)'^o,o = 0, then the function vri(x), as x approaches to its dominant singularity, has one of the three 
types of asymptotic properties shown in Case 1 to Case 3 of Theorem \4.^ 

Proof. First suppose that h 2 {x*,YQ{x*))'K 2 {Yo{x*))+hQ{x*,YQ{x*))-KQfl = 0, but/ii(Xo(yo), yo)7ri(Xo(yo))+ 
/io(Xo(yo)) yo)7ro,o / 0. We then have the following four cases: 

1 . li x* < xi < X 3 , then xi is a pole and the dominant singular point of 7 ri(x) since x* is a removable 

singular point of tti{x), which leads to ~ T^iix) = C, the same type in Case 1. 

2. If X* <xi = x^, the same type of asymptotic result as in Case 2 can be obtained. 

3. If X* = X 3 < ail, then the factor a/1 — xjxdom is cancelled out from both the denominator and 
the numerator in the expression for ' 7 ri(x). By considering 7 r^(x), we obtain the same type of asymptotic 
result as that given in Case 3. 

4. If x* = xi, then x* would be a pole of tt 2 {Yq{x)). This would imply 712 (^ 0 ( 3 :*)) = 00 , which 
contradict to h 2 {x*,YQ{x*))Tr 2 {Yo{x*)) + ho{x*,Yo{x*))7ro^o = 0 since h 2 {Xo{Yo{x*),YQ{x*)) = 0 implies 
h 2 {x* ,Yq{x*)) 7 ^ 0. Hence this case is impossible. 

Next, assume that both h 2 {x*,Yo{x*))7r2{Yo{x*))+ho{x*,Yo{x*))7ro^o = 0 and/ii(Xo(yo), yo)7ri(Xo(yo))+ 
/io(Xo(yo)) yo)'^o,o = 0- We then have the following two cases: 

1. If max{a:i, x*} < X 3 , then both aii and x* are removable poles of 7ri(x). By considering 7r^(x), 
we obtain the same type of asymptotic result as that given in Case 3. 

2. If maxjxi, x*} = X3, then the factor -^1 — xjxdom is cancelled out from both the denominator 
and the numerator in the expression for ti 2 {Yq{x)) if xi = X3 and for 7ri(x) if x* = X3. By considering 
7r^(x), we obtain the same type of asymptotic result as that given in Case 3. 

Finally, the case in which hi{Xo{yo),yo)'Ki{Xo{yo))+ho{Xo{yo),yo)'Kofl = 0, but h 2 {x*,Y o{x*))tt 2 (Yo{x*)) 
+/io(x*, yb(37*))vro,o / 0 can be similarly considered. □ 
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Remark 4.9 We believe that both h 2 {x*,YQ{x*))Tr 2 {YQ{x*)) + hQ{x*,YQ{x*))TTQfi / 0 and /ii(^o(yo)) yo) 
7ri(Xo(yo))+ /io(-^o(yo)) yo)'^o,o / 0 always hold, though at this moment we could not find a proof. How¬ 
ever, no new type of asymptotic property will appear without this condition as shown in Corollary If. ^ 
In the rest of the paper, the analysis will be carried out with this condition, which is also valid without 
this condition. 

Remark 4.10 When x^om = < min{x*, 5 ;i}, the numerator in the expression for tt\{x) is not zero 

at X 3 . 


5 Tail Asymptotics of Boundary Probabilities tt^ 0 and tto^u 

Since iti{x) and '/r 2 (y) are symmetric, properties for 7 ri(x) can be easily translated to the counterpart 
properties for 'n 2 {v)- Therefore, tail asymptotics for the boundary probabilities 7ro,n can be directly 
obtained by symmetry. 

The exact tail asymptotics of the boundary probabilities vr^^o is a direct consequence of Theorem 14.81 
and a Tauberian-like theorem applied to the function t^i[x). Specifically, if t^i{x) has only one dominant 
singularity, then Theorem 14.11 is applied; and if iifix) has two dominant singularity, then Theorem 14.21 
is applied. 

The following theorem shows that there are four types of exact tail asymptotics, for large n, together 
with a possible periodic property if vri [x) has two dominant singularities that have the same asymptotic 
property. 

In the theorem, let Xdom be the positive dominant singularity of 'ni{x). Consider the following four 
cases regarding which of x* , xi and X 3 will be Xdom '■ 

Case 1. Xdom = min{x*,xi} < X 3 with x* 7 ^ xi, or Xdom = xi = x* = X 3 ; 

Case 2. Xdom = X 3 = min{x*,xi} with x* / xi; 

Case 3. X 3 = Xdom < min{x*,xi}; 

Case 4. Xdom = x* = xi < X3. 


Theorem 5.1 Consider the stable non-singular genus 1 random walk. Corresponding to the above four 
cases, we have the following tail asymptotic properties for the boundary probabilities 'Kn,o for large n. In 
all cases, CQ^fixdom) A < * < are given in Theorem \4.^ 


1 . If Pi j is not X-shaped, then there are four types of exact tail asymptotics: 
Case 1: (Exact geometric decay) 

/ I ^ 

^n,0 ^ CQ^\{xdom) 


Xdom 


(5.1) 


Case 2: (Geometric decay multiplied by a factor of n ^l‘^) 




C0,2{xdom) ^-1/2 f _ \ 


IT 


Xdom 


n—1 


(5.2) 
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Case 3: (Geometric decay multiplied by a factor of n 


7rn,0 ~ 


Co,3 (^^dom) ^—3/2 ( ^ 


vr 


^dom 


n—1 


Case (Geometric decay multiplied by a factor of n) 


^n,0 ^ CO,4(^dom)^ 


^dom 


n—1 


(5.3) 


(5.4) 


2. If pij is X-shaped, but both p\^j and p\^j are not X-shaped, we then have the following exact tail 
asymptotic properties: 

Case 1: (Exact geometric decay) It is given by h5.1\) : 

Case 2: (Geometric decay multiplied by a factor ofn~^^‘^) It is given by 
Case 3: (Geometric decay multiplied by a factor of 


TTn.O ~ 


IfOjSi^dom) 4“ ( 1) Co,3( ®dom)] —3/2 


TT 


^dom 


n—1 


(5.5) 


Case 4- (Geometric decay multiplied by a factor of n) It is given by 

3. If pij and p[^j are X-shaped, but pf'^- is not, we then have the following exact tail asymptotic prop¬ 
erties: 

Case 1: (Exact geometric decay) When x* >xi, it is given by 115.1\} : when Xdom = x* <xi, it is 
given by 

/ ^ Nn-l 

^n,0 ^ [co,l (®dom) 4“ ( 1) Co,l( | ) , (^’S) 

V ^dom ) 


Case 2: (Geometric decay multiplied by a factor of n ^1“^) When x* > xi, it is given by ^5.21) : 
when Xdom = x* < xi, it is given by 


T^nfl ~ 


[co,2(3Jdom) 4“ ( 1) Co,2( —1/2 


TT 


^dom 


n—1 


(5.7) 


Case 3: (Geometric decay multiplied by a factor of n ^/^) It is given by US. 5)) . 

Case 4-' (Geometric decay multiplied by a factor of n) It is given by S- 

4. If pij and p^^j are X-shaped, but is not, then it is the symmetric case to 3. All expression in 3 

are valid after switching x* and xi. 

5. If all Pi j, p^^] and pf] are X-shaped, we then have the following exact tail asymptotic properties: 

Case 1: (Exact geometric decay) When x* <xi, it is given by 115.6\} : when x* > xi, it is also 
given by L5. by replacing the dominant singularity x* by xi. 

Case 2: (Geometric decay multiplied by a factor of When x* < xi, it is given by (B; 

when X* > xi, it is also given by (B by replacing the dominant singularity x* by xi. 

Case 3: (Geometric decay multiplied by a factor of n~^/^) It is given by 115. 51) . 
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Case 4- (Geometric decay multiplied by a factor of n) It is given by 

/ 1 

~ [co,i{Xdom) + (-l)”“^Co,4(-a:dom)] fl I - ) . (5.8) 

\ ^dom J 

Proof. 1. Since Pi^j is not X-shaped, all —xs, —x* and —xi are not dominant singularities according 
to Corollary 13.11 Theorem 14.51 and Theorem 14.61 Therefore, there is only one dominant singularity for 
7 ri(x). The tail asymptotic properties of follow from Theorem 14.81 and the direct application of the 

Tauberian-like theorem (Theorem 14.ip . 

2. We only provide a proof to the cases, which are not identical to that in 1. 

Case 1. For the case that x^om = = x* = 3 : 3 , we notice that —X 3 is also a dominant singularity 

1 Corollary 13.ip . In this case, the Tauberian-like theorem (Theorem 14.21) is used to have a tail 
asymptotic expression consisting of two terms, one, corresponding to the positive dominant singu¬ 
larity, with the exact geometric decay rate and the other, corresponding to the negative dominant 
singularity, with the geometric decay rate multiplied by a factor of Therefore, the term with 

the geometric decay rate is the dominant (decay slower) term leading to the same tail asymptotic 
property given in (j5.1l) . 

Case 2. Similar to Case 1, —X 3 is also a dominant singularity. The Tauberian-like theorem (Theo¬ 
rem leads to a tail asymptotic expression consisting of two terms, one with the geometric rate 
multiplied by a factor of (dominant term) and the other by 

Case 3. In this case, both X 3 and —X 3 are dominant singularities having the same asymptotic property 
according to Theorem 14.81 The tail asymptotic expression follows from the application of the 
Tauberian-like theorem (Theorem 14.2p . 

3. In this case, —X 3 and —x* are singularities, but —xi is not. We only provide a proof to the cases, 
which are not identical to that in 1 or in 2 . 

Case 1. For the case when x* = xi = X 3 , there are two dominant singularities. The Tauberian-like 
theorem (Theorem I4.2jl leads to a tail asymptotic expression consisting of two terms, one (corre¬ 
sponding to the positive singularity) with a geometric decay rate, and the other (corresponding 
to the negative singularity) with the same geometric decay rate multiplied by a factor of 
that is dominated by the geometric decay. 

When X* < xi, both x* and —x* are dominant singularities with the same asymptotic property, 
which leads to the tail asymptotic expression by using Theorem 14.21 

Case 2. For case when X 3 = x*, there are two dominant singularities having the same asymptotic 
property. The tail asymptotic expression follows from Theorem 14.21 

Case 4. In this case, there are two dominant singularities, but the contribution from the positive 
dominant singularity dominates that from the negative dominant singularity. The tail asymptotic 
expression follows from Theorem 14.21 

4. The symmetric case to 3. 

5. In this case, all —x*, —x* and —X 3 are singularities. We only provide a proof to the cases, which 
are not considered in the above. 
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Case 1. The only new situation here is the case when x* = x* = X 3 . In this case, we have the same 
asymptotic property at both dominant singularities, which leads to (15.6p . 

Case 4. In this case, we have the same asymptotic property at both dominant singularities, which 
leads to (15.81) . 


□ 

From the above theorem, it is clear that if there is only one dominant singularity, then the boundary 
probabilities tt^^o have the following four types of astymptotics: 1. exact geometric; 2. geometric 
multiplied by a factor of 3. geometric multiplied by a factor of and 4. geometric 

multiplied by a factor of n. If there are two dominant singularities, but with different asymptotic 
properties, 7r„ q also has one of the above four types of tail asymptotic properties. Finally, if we have 
the same asymptotic property at both dominant singularities, then vr^^o reveals a periodic property with 
the above four types of tail asymptotics, which is a new discovery. 


6 Tail Asymptotics of the Marginal Distributions 


In the previous section, we have seen that the asymptotic behaviour of the function 7 ri(x) ( 7 r 2 (y)) 
at its dominant singularity or singularities determines the tail asymptotic property of the boundary 
probabilities tt^^o ('^o,n)- According the the fundamental form of the random walk, it, together with the 
property of the kernel function h{x,y), also determines the tail asymptotic property of the marginal 
distribution TTnj (and = Yli 

In this section, we provide details for the exact tail asymptotics of the marginal distribution 

( 2 ) 

The exact tail asymptotics of ^ can be easily obtained by symmetry. First, based on the fundamental 
form, we have 

hi{x,y)TTi{x) + h2{x,y)TT2{y) + ho{x,y)TTofi 

TT =-tt-^- 

-h{x,y) 

and therefore, 


7 r(x, 1) 


hi{x, l) 7 ri(a;) + h 2 {x, l) 7 r 2 (l) + hpjx, I) 7 ro,o 
—h{x, 1 ) 

hi{x, l) 7 ri(x) + h 2 {x, l)vr2(l) + ho{x, l)7ro,o 

-S(l)[x-Xo(l)][x-Xi(l)] 


If Mx > 0, then Aii(l) = I, which implies that the denominator of the expression for 7 r(x, 1) does not 
have any zero outside the unit circle. In this case, has the same tail asymptotics as vr^^o- The only 
difference is the expression for the coefficient, which can be obtained from straight forward calculations. 

If Mx < 0, then Ao(l) = I and Ai(l) > I. Ifpi,j is not X-shaped, the analysis is so-called standard, 

(k) 

details of which will be provided here. If pij is X-shaped, then there are four subcases based on if p) ■ 
is X-shaped or not. For these cases, detailed analysis varies, but similar. We provide details here for 
the case where both for A: = 1, 2 are not X-shaped. Let 2 ; = min{x*, xi}. and consider the following 
four cases: 

1. min(Xi(I),z) < X 3 and Ai(l) / z. In this case, has an exact geometric decay with the 
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decay rate equal to Xdom = min(Xi(l), z): 


r(i) ^ 

V Xdom 


n—1 


where 


(x) _ 


Cl = 


[/il(Js:i(l),l)7ri(Xi(l))+/i2(Xi(l),l)7r2(l)+feo(Xi(l),l)7ro,o]Xi(l) ^ ^ ^ 

-S(l)(Xi(l)-l)-> < Z, 


hx{z,l)co^^{z)+h2{z,l)TT2(l)+ho{z,l)-KQfi 
-H(l)h-Xo(l)]h-Xi(l)] 


^l(l) > 


with Co, 1 ( 2 ;) being given in Theorem! 

2. Xi(l) = 2 ; < X 3 . In this case, Xi(l) = Ti is impossible, since otherwise /i(Xi(l),l) = 0, which 
implies 1 = lo(xi) or 1 = Yi(xi). This is contradiction to To(51i) > 1- Hence, only Xi(l) = x* may 
hold. There are two subcases: 

2(a): 1 = In this case, hi{x*,l) = hi{x* ,Yq{x*)) = 0. We can write /ii(x, 1) as /ii(x, 1) = 

ai{x) + bi{x) = (x — Xi(l))/i^(x) with being a linear function of x, which yields 


7r(x, 1) = 


hi{x, l)7ri(x) + h 2 {x, l)7r2(l) + hpjx, l)7ro,o 
-a(l)(x - l)[x - Xi(l)] 
hl{x)Tri{x) ^ h2{x, l)7r2(l) + ho{x, l)7ro,o 


—a(l)(x — 1) —a(l)(a: — l)[x — Xi(l)] 

Therefore, 7r{x, 1) has a single pole Xi(l), which leads to an exact geometric decay (recalling 7 r(l, 1) 7 ^ 1): 


( 1 ) ^ ^ 
Xi{l) 


'^n ~ ^2,1 


n—1 


with the coefficient given by 


41= lim 1-^ 


7r(x,l) = 


[h2(^i(l), l)7r2(l) + ho(Xi(l), l)7ro,o]Xi(I) 
h(l)(Xi(l)-l) 


2 (b): 1 = Yi{x*)- In this case, hiix* ,Yq{x*)) = 0 and 10(2:*) < Yi(x*), we obtain hi{x*,l) = 
hi{x*,Yi{x*)) > 0, which implies that x* is a double pole of vr(x, 1) (noting that /i 2 (x*, l) 7 r 2 (l) + 
ho(x*, l)7ro,o > 0 since h 2 (x*,l) > /i 2 (Xo(l),l) = 0 and /io(x*, l)7ro,o > 0). The corresponding tail 
asymptotic is given by 

(1) (*) f ^ 

where 


n—1 


c(^) = lim I 1_ — 

2.2 Xi(I) 


7r(x, 1) 


X-J-Xi(l) \ 

Xi(l)[hi(Xi(l), l)co,i(x*) + h 2 (Xi(l), l)7r2(l) + ho(Xi(l), l)7ro,o] 

S(I)[Xi(l)-l] 


with co,i(x*) given in Theorem W. 

3. min(Xi(l), z) = X 3 . In this case, there are four possible subcases, for which proofs are omitted 
since they are similar to that for the previous cases: 
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3(a): Xi(l) = z = X 3 leading to an exact geometric decay: 




n—1 


where 


cl*? = lim (1 -^ = 


/i2(x3, l)vr2(l) + ho{x3, l)7ro,o 


X3a{l){x3 - 1 ) 

3(b): Xi(l) = X 3 < z leading to an exact geometric decay: 


n—1 


where 


7r(i) ~ ( — 


^(x) ^ A _ ^ hi{x 3 , l)7r(x 3) + h 2 {x 3 , l)7r2(l) + ho{x 3 , l)7ro,o 

X^X 3 I X 3 J ’ 


X3a{l){x3 - 1) 

3(c): z = X 3 < -^i(l) with x* / aii leading to a geometric decay multiplied by the factor 


(l") (x) - 1/2 


X3 


n—1 


where 


(x) 

Co q = lim 

X— 1 x 3 


f. / 11 hi(x3,l)co,2(x3} 


S(l)(Xi(l)-l)[Xi(l)-X3] 

with Co ,2 ( 2 : 3 ) given in Theorem 14.81 

3(d): z = X* = xi = X 3 < -Ti(l) leading to an exact geometric decay 


r(i) ^ (1 


'^n ~ C 3 4 


X3 


n—1 


where 




hi{z, l)co,i(z) + h 2 {z, 1 ) 712 ( 1 ) + ho{z, 1 ) 710,0 

-S(l)[z-Xo(l)][z-Xi(l)] 


4. X 3 < min(z,Xi(l)) leading to a geometric decay multiplied by the factor n 


71, 


( 1 ) 


~ C, 


(*)^-3/2 f ]_ 


X3 


n—1 


where 

C 4 *^ = lim (\ —7r'(x, 1) 
X^X 3 \ X 3 J 

with Co, 3 (^ 3 ) given in Theorem 14.81 


hi { x 3 , 1)00,3(3:3) 
a(l)(a:3 - l)[-’^i(l) - X 3 ] 


For the completeness, we provide a summary of tail asymptotic properties for the marginal distri¬ 
bution Tii^^ for all possible cases. For this purpose, let Xdom be the positive dominant singularity of 
7 r(x, 1). Note that Xi(l) 7 ^ xi. The following are the all possible cases according to which of xi, x*, X 3 
and Xi(l) is Xdom- 
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Case A. Xdom = min{xi,x*,X 3 } < Xi(l); 
Case B. Xdom = < min{xi,x*,X 3 }; 

Case C. Xdom = Ai(l) = x* < min{xi,X 3 }; 
Case D. Xdom = Xi{l) = X 3 < x*\ 

Case E. Xdom = =X3 = x*. 


Remark 6.1 The cases here are different from the cases classified in the previous section and the next 

section. 

The exact tail asymptotic properties are obtained according to the expression of 7r(x, 1) and the 

Taubarian-like theorem. 

Theorem 6.1 For the stable non-singular genus 1 random walk, the exact tail asymptotic properties 

for the marginal distribution 7rn \ as n is large, are summarized as: 

Case A: This case includes Cases I -4 in the previous section. iTn^ has the same types of asymptotic 
properties as vr^^o given in Theorem 15.11 respectively, with possible different expressions for the 
coefficients. 

Case B: TTn'^ has an exact geometric decay. 

Case C: has an exact geometric decay ifYo{x*) = 1 and a geometric decay multiplied by a factor 

of n ifYi{x*) = 1 , respectively. 

Case D: TTn'^ has an exact geometric decay. 

Case E: TTn'^ has an exact geometric decay. 


7 Tail Asymptotics for Joint Probabilities 


In the previous sections, we have seen how we can derive exact tail asymptotic properties for the 
boundary probabilities and for the marginal distributions based on the asymptotic property of '7ri(x) 
(^ 2 ( 2 /)) and the kernel function. However, the exact tail asymptotic behaviour for joint probabilities 
cannot be obtained directly from them. Further tools are needed for this purpose. Our goal is to 
characterize the exact tail asymptotics for for each fixed j and for each fixed i. Due to the 
symmetry, in this section, we provide details only for the former. 

The relevant balance equations of the random walk are given by 


(1 -Po?o)^o,o 
(1 -p[)!o)^i.o 
(1 - PoIo)’ri.o 
(1 “ 


( 1 ) , ( 2 ) , 

Pli, 0 ^ 1 ,0 +Po,- 1 ^ 0,1 +P- 1 -iTTgl, 

( 0 ) I ( 1 ) I , ( 2 ) , 

Pi,oTTo,0 +Pli,0^2,0 + P-1,-I7r2,l +Pi,li7ro,l +P0,-l7ri,l, 

+ P-1,-I7ri+I,i + Pi,_i7ri_i,i + po,-i7ri.i, i > 2 , 
+P-l,-l^i+l,j+l +P0,-l^i,j+l + PlfiT^i-lJ + P-lfiT^i+lJ 
+Pl.lT^i-l,j-l + P0,l'^i,j-1 + j > 2 . 
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Let 


i=l 

oo 

V’i(y) = '^'^i,jy''~^, * > 0. 
i=i 

From the above definition, it is clear that ipo{x) = tti{x) and V'o(y) = ^ 2 ( 2 /)- From the relevant balance 
equations, we obtain 


or 


where 


c{x)^pi{x) + hi{x)ipo{x) = al{x), 
c{x)ip 2 {x) + b{x)ipi{x) + ai{x)ipo{x) = a*(x), 
c{x)ipj+i{x)+ h{x)Lpj{x) + a{x)Lpj-i{x) = a*{x), j>2, 


(fj+i{x) = 


—b{x)ipj{x) — a{x)ipj-i{x) + a*{x) 
c(x) ' 


j >0, 


ao(®) = -C2(a;)vro,i - bo(x)7ro,o, 

a’^(x) = -C2(x)7ro,2 - f'2(a:)vro,i - ao(x)7ro,o, 

a*(x) = -C2(x)7roj+i - b2(x)7roj - a2(x)7roj-i, j > 2. 


(7.1) 

(7.2) 

(7.3) 


(7.4) 


First, we establish the fact that a zero of c{x) is not a pole of ipj{x) for all j > 0. Therefore ^j{x) 
has the same singularities as ^Pq[x). 

Let y = Yq{x) be in the cut plane C^,, and let ydom and Xdom be the positive dominating singular 
points of ijjoiy) and ^po{x), respectively. Let 


OO 0000 

fk{x) = -a 2 {x) - & 2 (x) - C 2 (x) k>l, 

j=k-l j=k j=k+l 

then, 


fi{x) = yf 2 {x) - C2(x)7ro,2 - b2{x)7ro^i, 

fk{x) = yfk+i{x) + a*k{x), k>2. (7.5) 

According to Theorem 14.31 when |x| < Xdom, we obtain 

hi{x,y)ipoix) = -/i2(x,y)V’o(2/) - ho{x,y)TTofi 

= y[fi{x) - ao(x)7ro,o] + a^ix) (7.6) 

= y^ f 2 {x) + yal{x) + a^ix) (7.7) 

= y^h{x) + y‘^al{x) + yal{x) + ao(x). (7.8) 


Let u = Since a zero of c{x) is a zero of Yq{x), u is analytic on the cut plane Cx- Using 

^ = -yi(x) - Yo{x) and ^ = yi(x)yo(a;), we obtain 

1 + b{x)u = —yua{x) and —= —a{x)u^. (7.9) 
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Write a = a(x), h = 6(x), c = c(x), a* = ai{x), bi = 6j(x), a* = a*(x), fj = fj{x) and ipj = ^j{x). 
We have following Lemma, which confirms that a zero of c(x) is not a pole of <fj{x) for all j > 0. 
Therefore, ipj{x) has the same singularities as ipo{x). 


Lemma 7.1 Let 

W-i = -, Wo = —bi, Wj = buwj-i + (1 + bu)wj-2- 

au 

Then, 

(-1)-^ [buwj-i + (1 + bu)wj- 2 ] + bi+ aiy = (-l)-^(l + bu)wj-i, j > 1, 

and 

hupi = yuf2Wo + ugi, 

j-2 

hiifj = {-iy^^yufj+iWj-i + u^{-iy^^~^a*_i,Wj_i_k{au)^ + ugi{auy~\ 

k=0 

where gi = aQ(x)oi(x) — hi{x)a\{x). 


(7.10) 

(7.11) 


(7.12) 
j>2, (7.13) 


Proof. By applying (|7.9I) and (I7.10p . we easily obtain equation (17.lip for j = 1. Assume that equation 
(17.lip is true for j < k, we show 

(_l)fc+i + (1 + bu)wk-i] +bi + aiy = (-1)*'+^(1 + bu)wk. (7.14) 

From the inductive assumption and the definition of wj, we have 

bi + aiy = + bu)wk-i[buwk-i +{I+ bu)wk-2] , 

= {-1^(1+ bu)wk-i +i-l)^^^Wk, (7.15) 

which yields equation (I7.14p . Equation (I7.12p is obtained by the direct substitutions of equations (|7.ip 
and (17.61) . 

Next, we show equation (j7.13p . We use the induction again. According to equations (j7.2h . (j7.7p and 

dnn), 


c{x)hiip2 = —bhupi — aihiifo + hial 

= -b[yuf 2Wo + ugi] - ai[y^f2 + yal + Oq] + al[aiy + bi]. 


It follows from equations (|7.9I) and (|7.10l) that 


hnp2 = 


which gives equation 


-bu^f 2 Wo - aiyuf 2 + u{au)gi 


-uf2 


' l + bu 

buwQ + ai - 


+ u{au)gi 


—au 

—uf 2 [buwo + (1 + bu)w-i] + u{au)gi = —yuf^wi — a^uwi + u{au)gi, 

(I7.13P for j = 2. Assume that equation (I7.13p is true for j < n. We prove the 
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result for j = n + 1. From equations (I7.9p . (|7.11l) . (I7.5p and the inductive assumption, we have 
c{x)hiipn+i = -bhiifn - ahupn-i + hial 

n—2 

= {-l)'^^‘^ybufn+iWn-i + - bugi{au)^~^ 


k=0 

n —3 


+{-l)^+^yaufnWn-2 + ou^(-l)"+^ ^al_i_kWn-2-k{au)^ - gi{au)'^ ^ + a;^[aiy + 6 i] 


k=0 


= {-l)'^~^‘^yfn+i[buwn-i - yauWn-2] + < {{-lT~^‘^buWn-i + yauWn-2 + aiy + 61 } 


n—2 


n—3 


+bu ^(-1)”+^ ’"a*^_j^Wn-i-k{auf + ^(-1)”+^ - (1 + bu)gi{au) 


n—1 


k=l fc=0 

= {-lY^‘^yfn+i[buWn-i + (1 + bu)wn- 2 \ + (-l)"+^a* (1 + bu)Wn-l 

n—2 

\n—l 


k=l 


+(1 + bu) ^(-1)'"+^ ^a*^_kWn-i-k{au)^ - (1 + bu)gi{auY 
which yields 

hupn+i = {-lT^‘^ufn+iWn + a*^au^Wn-i + avY a*^_^Wn-i-k{auY' + ugi{auy 

k=l 

n—1 

= {-lY^‘^yufn+iWn + U + ugiiau)"^. 

k=0 

This completes the proof. □ 


n—2 


Corollary 7.1 ifo{x) and (pj{x), j > 1, have the same singularities. 

The following Lemma is useful in characterizing the tail asymptotics of vr^j- for a hxed j. 

Lemma 7.2 //min{x*,xi} > X 3 , then 

lim . / 1 - — (f', (x) = C 3 J {Xdom ), 

y Xdom 

where c^^o{xdom) ^5 given in Theorem \4.8\ and 

(^dom) — [^si^dom) “1“ ^3(^dom)j] ( ; J ^ 0, 


with 




i^dom) 


^3,0(^dom)^! i^dom'} 

^i^dom) 

hi {Xdom^ ^(^dom))^3,o(^dom) 


c(^dom) 


(7.16) 

(7.17) 

(7.18) 
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Proof. When min{x*, Xi} > X 3 , we have Xdom = ± 2 : 3 . Without lose of generality, we assume Xdom = 
in the proof. Since ^j{x), j > 0, is continuous at X 3 , — -^^jix) = 0. Let j = 1. Then, 

/, ^ -c’{x)ipi{x) - bi{x)(p'o{x) - bi{x)ip'o{x) + a^'{x) 

=-Jw- 

and 


V ^ ^ -&i(a;3)c3,o(a;3) , ^ 

hm . I -= - - -r- = C 3 ,l(X 3 ). 

x^X3 y X 3 c(X 3 ) 


Assume that Imix^xz y^l — ^^'k{x) exists for k < j and 

/ X 

lim Jl- (p'f.{x) = C 3 ,fc(x 3 ), 

x^x3 y X3 


we obtain 


(f'k+iix) = 


-c'{x)ipk+i{x) - b{x)ip'f.{x) - b'{x)ipk{x) - a{x)ip'f._^{x) - a'{x)ipk-i{x) + a*f!{x) 


and 


lim . /1 —(x) = 

x^xay X 3 ^ 


X—>X^ 

Therefore, we can inductively have 


c{x) 


-b{x3)c3^k{x3) - a{x3)c3^k-l{x3) 

C{X3) 


= C3,fc+l(x3). 


C3,l(a;3)c(x3) + C3,o(x3)6i(x3) = 0, 

C3,2(a;3)c(a:3) + C3,i(x3)6(x3) + 03,0(2:3)01 (X3) = 0, 

C3j+i(a:3)c(x3) + b{x3)c3j{x3) + a(x3)c3,j_i(x3) = 0, j > 2. 


(7.19) 

(7.20) 

(7.21) 


It follows that {c 3 ,fc(x 3 )} is the solution of the second order recursive relation determined by equations 
(j7.19p - (l7.2ip . Since b‘^{x 3 ) — 4 a(x 3 )c(x 3 ) = 0, C 3 j{x 3 ) takes the form given by equation (17.161) . ^ 3 (^ 3 ) 
and i? 3 (x 3 ) are obtained by using the initial equations: 


[^ 3 ( 2 : 3 ) + B3{x3))c{x3)] 
Yi{x3) 


A3(x3)c(x3) +C3,o(x3)6i(x3) = 0, 

+ A3(x3)6(x 3) + C3,o(x3)ai(x3) = 0. 


We are now ready to prove the main theorem of this section, in which 

Cl o(^c/om)^l i^dom) 




^i^dom) 

^1 {^domj ^(^dom)) 


□ 


(7.22) 


^l{^dom}\ ( \ 

7 ^ I ^1,0 \^dom) ^ 


C'(^dom)[^1 (^dom) ^^(^dom)]^(^dom) ^i^dom) J 

C2,o(^dom)&l i^dom) 


^si^dom) is given in (|7.17p . 


'^2(^dom) — 

"^4(^dom) 


c(3^dom) 

(^dom)co,4 i^dom) 
c(3^dom) 


(7.23) 

(7.24) 
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B\{X(iom) — 


'®2(^dom) 

and B 3 {xdom) is given in (|7.18p . 


hi (Xdomi ^(®dom))ci,o(^(iom) 

1^1 i^dom') ^oi^dom)]^oi^dom') 
(^2fl{Xdom')hl{Xdomi ^oi^dom)) 


OiYq {Xdonri)^ 


(7.25) 

(7.26) 


Theorem 7.1 Consider the stable non-singular genus 1 random walk. Corresponding to the four case, 
we then have the following tail asymptotic properties for the joint probabilities Tin,] for large n. 

1. If Pi j is not X-shaped, then there are four types of exact tail asymptotics: 

Case 1: (Exact geometric decay) 


iXdom) 


1 


(p^dom) 


i-1 


-|- Bi (xdom) 


1 


^0 (^dom) 


J-I 


1 


^dom 


n—1 


Case 2: (Ceometric decay multiplied by a factor of n J 

[A2ixdom) + {j - l)B2{xdom)] f 1 -1/2 f ^ ^ 

^ 0F \Yi{xdom)) "" \xdom) 

Case 3: (Ceometric decay multiplied by a factor of n~‘^l’^) 


j > 1; (7.27) 


n—1 


j > 1; 


(7.28) 


[A^(^Xdom) “1“ {j ^')Bn^Xdom)] 


TT 


. Yl {xdom) 

Case 4' (Ceometric decay multiplied by a factor of n) 


1 \/I 
1 ^ n-3/2^ ^ 


Xdom 


n—1 


j > 1; 


A4 {Xdom ) 


1 


{xdom') 


3-1 


n 


Xdom 


n—1 


, j>l- 


(7.29) 


(7.30) 


2. If Pij is X-shaped, but both p\^j and pf'j are not X-shaped, we then have the following exact tail 
asymptotic properties: 

Case 1: (Exact geometric decay) It is given by 

Case 2: (Geometric decay multiplied by a factor of n~^/^) It is given by Ii7.28\ j; 

Case 3: (Geometric decay multiplied by a factor of n~‘^l‘^) It is given by 


[AsiXdom) + Asi-Xdom)] 


TT 


1 


Yl {xdom} 


i-1 


n 


-3/2 


^dom 


n—1 


, j > 1; 


(7.31) 


Case 4' (Geometric decay multiplied by a factor of n) It is given by \7.30{) . 


38 
































3. If pij and p\^j are X-shaped, but pf'^- is not, we then have the following exact tail asymptotic 
properties: 

Case 1: (Exact geometric decay) When xi < x*, it is given by {1.2% ; when xi = x* = x^, it is 
also given by l7.2''/\) ; when x* <xi, it is given by 


'^n,j ~ [^l(3Jdom) “1“ ( 1) '^^l( 


1 


i-1 / 1 ^ 


J \Xdom 


, i > 1; 


(7.32) 


Case 2: (Geometric decay multiplied by a factor of n When x* > xi, it is given by Ii7.28( ); 
when X* <xi, it is given by 


[A2{xdom) + A2{-Xdom)] f 1 V ^ ^-i /2 f 1 


TT 


^1 i^dom') 


n 


Xdom 


n—1 


, j>i; 


(7.33) 


Case 3: (Geometric decay multiplied by a factor of n J It is given by 
[Az{Xdom) + A^{-Xdom)] { 1 -Z/2 ( ^ 


n—1 


\^li^dom) J \^dom 

Case 4-' (Geometric decay multiplied by a factor of n) It is given by 17.30{) . 


, j > 1; 


4- If Pi j and p^^] are X-shaped, but p^P- is not, then it is the symmetric case to 3. All expression 
in 3 are valid after switching x* and xi. 

5. If all Pi i, p^^] and pf'] are X-shaped, we then have the following exact tail asymptotic properties: 

Case 1: (Exact geometric decay) When x* <xi, it is given by 1 7. 32\ ); when x* > xi, it is also 
given by (f 7 . 320 replacing the dominant singularity x* by xi; 

Case 2: (Geometric decay multiplied by a factor of n~^l’^) When x* < xi, it is given by (f 7. 33^ ; 
when X* > xi, it is also given by 1 7. 33t) by replacing the dominant singularity x* by xi; 

Case 3: (Geometric decay multiplied by a factor of n~^^'^) It is given by \7.31^ ; 

Case 4' (Geometric decay multiplied by a factor of n) It is given by 


'Xn,j ~ [^4(3^(iom) “1“ ( 1) ^A/^(^ ^^ciom)] 


1 


^1 {Xdom) 


i-1 


n 


^dom 


n—1 


j > 1- 


Proof. 1. 

Case 1: It follows from Section 021 that { 1 — ) 7^o{x) = cq i(xdom)- By the induction 

aom Y ^dom J ^ ^ ' 

and equations ([10])-([731), lirnx^xdom = ci,j{xdom) with 


Cl,1 {xdom)dom) P g(^dom)^l(®dom) 

Cl,2iXdom)c{Xdorn') P Ci^ii^Xdom^bi^Xdom) P Cl^o(®dom)®l(^dom) 
ClJ-\-l{xdom)dxdom) P Ci^j{xdom)b{xdomt) P Cij—i(^Xdom)rL{xdom) 


0 , 

0 , 

0, j>2. 
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Since cij{xdom)^ 3 ^ 0, satisfies the second order recursive relation above, it takes the form of 




{.^dom) 


-|- B\ {Xdorri) 


^ iS^dom) 


j >0. 


To determine Ai = Ai{xdom) and Bi = Bi{xdom)^ use the initial equations: 


(^1 T T ^1,0(^dom)^l (^dom) 0? (7.34) 

^ c{^dom) T (-^1 T Bi^bi^Xdom') T Cl,o{^do7n')^l{^doni^ — 0* (7.35) 

Multiplying both sides of equation f|7.35p by Yo{xdom)^ adding the resulting one to f|7.34p . and taking into 

account CL{Xdom}YQ {Xdom) T ^{^dorn)Y()(^Xdorn} T ^(^dom) 0; (^dom; (^dom)) ^l{^dorn}Y()(^Xdorn) T 

bi{xdom^ and ci^Xdom) — Y()(^Xdom)Yi(^Xdom)^{^dom') yield. 




Yl {Xdom) 




Yq {Xdom) 


(■^1 Y Bi)c(^Xdom) T ^1,0(^dom)^l (^dom) 0; 


Yq {XdoTTi) 


Y (x ^^(^dom) T Bic(^Xdom) T (-^1 Y Bi^b(^Xdo7n)YQ(^Xdom) Y Ci^Q(^Xdom)^l{^do7n)Y{)(^Xdo7n) 0; 

which gives (|7.25p and (|7.22p . So, Bi{xdom) = 0 if Xdom = x* and Bi{xdom) / 0 if Xdom = By the 
Tauberian-like theorem, we obtain f|7.27p . 

Case 2: Similar to that for 1-Case 1. From the proof, we have (|7.26p and (|7.23p . 

Case 3: Write 

oo oo / \ n 

+ l)7rn+2,jX” = + l)x^7rn+2,j { — ) ■ 

n=0 —n 


n=0 


According Lemma 17.21 and the Tauberian-like theorem, we have 

/ , n C3j(a;3) _i/2 

[n + l}x^7rn+2,j ~ ^ n ' , 


which is equivalent to (I7.29p . 

Case 4: The results can be proved in the same fashion as in Case 1 and Case 2. 

The proofs of the other cases are omitted due to the similarity to 1 and Theorem 15.11 □ 


8 Examples and Concluding Remarks 

In this paper, for a non-singular genus 1 random walk, we proposed a kernel method to study the exact 
tail asymptotic behaviour of the joint stationary probabilities along a coordinate direction, when the 
value of the other coordinate is fixed, and also the exact tail asymptotic behaviour for the two marginal 
distributions. A total of four different types of exact tail asymptotics exists. The fourth one, a geometric 
decay multiplied by a factor n, was not reported before for this discrete-time model (the same type was 
reported recently for a continuous-time random walk model by Dai and Miyazawa my In this study, 
we also revealed a new periodic phenomena for all four types of exact tail asymptotics when there are 
two dominant singularities for the unknown generating function, say 7ri(a;), with the same asymptotic 
property at them. 
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The key idea of this kernel method is simple and the use of the Tauberian-like theorem greatly 
simplifies the analysis, which, unlike in the situation when a standard Tauberian theorem is used, is 
also rigorous. Under the assumption that there is only one dominant singularity, this method provides 
a straightforward routine analysis for the exact tail asymptotic behaviour. However, without this 
assumption, the analysis is not simple, at least to our best effort, for telling how many dominant 
singularities and when a pole is simple. It is also challenging to characterize the exact tail asymptotic 
along a coordinate direction when the value of the other coordinate is not zero, since it is not a direct 
consequence of the kernel method. 

This kernel method can also be used for characterizing the exact tail asymptotics for the non-singular 
genus 0 case and the singular random walks (see Li, Tavakoli and Zhao [32]). With the detailed analysis 
provided in this paper, we expect further research in applying this kernel method to more general 
models. 

The complete characterization of the exact tail asymptotic behaviour provided in this paper does 
not necessarily imply that for any specific model, a characterization explicitly in terms of the system 
parameters exists. However, we are confident that for any specific model, if using a different method 
could lead to a such characterization, in terms of system parameters, then it can be done using the 
kernel method. Finally, we mention two examples, which have been analyzed by using the proposed 
kernel method. 

Example 1. A generalized two-demand model was considered in Li and Zhao [36| using the same 
idea proposed in this paper. For this model, let A and {k = 1,2) be the Poisson arrival rate with 
two demands and the arrival rate of the two dedicated Poisson arrivals, respectively. Furthermore, let 
Hk {k = 1,2) be the exponential service rates of the two independent parallel servers. For a detailed 
description of the model, one may refer to [36] • For this model, the three regions, on which the joint 
probabilities along a coordinate direction, say queue 1, have an exact geometric decay, a geometric 
decay multiplied by a factor and a geometric decay multiplied by a factor are extremely 

simple, which are: (a) ; (b) and (c) ^ respectively. 

(k) 

Example 2. Consider the simple random walk, or a random walk for which pij and both 
{k = 1,2) are cross-shaped. We then can follow the general results obtained in this paper to have 
refined properties. For example, consider the case of My > 0 and Mx < 0 and assume that the system 
is stable. Then, along the x-direction, has three types exact asymptotics in the following respective 
regions: 

1. Exact geometric: 


X3 


X3-I 


PO-1 

P0,1 


- 1 


( 1 ) , ( 1 ) ^ ( 1 ) 
Po,l +Pi,M3 >P-i,o; 


2. Geometric with a factor n 


X3 


X3-I 


Po-i 

Po,i 


- 1 


( 1 ) , ( 1 ) ( 1 ) 
Po,i +Pi>3 =p_i,o; 


3. Geometric with a factor n 


X3 


X3-I 


Po-i 

Po,i 


- 1 


( 1 ) , ( 1 ) / ( 1 ) 
Po,i +Pi,M 3 <P-io- 
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When My < 0 and < 0, this example also reveals the fourth type of exact tail asymptotic 
property, or a geometric decay multiplied by the factor n along the x-coordinate direction in the region 
defined by the following conditions: 


a:3 

X3 - 1 


Po-i 

Po,i 


- 1 


( 1 ) , ( 1 ) \ ( 1 ) 
P0,l +Pl>3 >P-l,o> 


( 8 . 1 ) 


and 


1/3 

1/3-1 


P- 1,0 

Pi,o 


- 1 


( 2 ) , ( 2 ) . ( 2 ) 
H,o+Po,i2/3 


hi{x*,yo) = 


P- 1,0 

Pi,o 


0 , 

Ji) 

P-i,o 

“oT’ 

Pifi 


( 8 . 2 ) 

(8.3) 

(8.4) 


{x* - l)p[,^liPi,o -1 _ ^ ^ {x* - l)[pL^j o 

(x* - l)Po}pi,o + pSpo,i P^a;* 


Here, x* € ( 1 , 3 ^ 3 ] and y* € (ljP 3 ] are the zero hi{x,YQ{x)) and /i 2 (Xo(y),y), respectively, whose 
existence is guaranteed by Lemma HlSl under conditions (18.ip and (18.2p : yo = 1 ^ 0 ( 3 :*) and in this case we 
have yo = y*] and xq = Xo(yb(3:*)). 


It is not very difficult to see this is not an empty region. The last thing which we need to check is 
the coefficient _ _ _ 

, X /i2(3;dom,2/*)[/ii(3;o,2/*)7r(xo) +/io(xo,p*)]vro,o] , „ 

CoAixdom) = - 


x*2 h[ [Xdom ,y*)yi [Xdom ) h'a (Xq (y*), y * ) 


( 8 . 6 ) 


or 


hi{xo,y*)'K{xo) + /io(xo,y*)7ro,o + 0, 
which is true since h 2 {xdom,y*) = h 2 {Xi{y*),y*) > h 2 {Xo{y*),y*) = 0. 
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A Proof to Lemmas I4.4H4.7I and Propositions I4.1H4.2 


Proof, of Lemma 14.41 1. From h{x,y) = 0, we have 

, _ a'{x)y^ + h’{x)y + c'{x) 

^ 2a{x)y + b{x) 

Using a(l) + 6(1)+ c(l) = 0, the property in (12.2p and the expression for in Lemma [3.11 we obtain 

(a) and (b). (c) is obvious. 

2. There are two possible cases: My < 0 and My > 0. If My < 0, according to the ergodicity 
condition in Theorem 12.11 M^'^M^ — MyM^'^ < 0 must hold, which yields 


/'(I) 


a(l)/i;(l,yo(l))/T(l,H(l)) 

a(l) [«;(!) + ai(l)yo'(l) + 6;(1)] /o(l,yi(l)) 
a(l)/ii(l,yi(l)) 


-M„ 




< 0 . 


From equation (14.8p . /(xs) > 0, it follows that f{x) = 0 has a root in (IjXs] since /(I) 

/'(I) < 0. 

If My > 0, we have 


/(I) = a(i)/ii(i,yo(i))/i;(i,yi(i)) 

-a(l)/ti(l,yo(l))[M^Mi^^ 

My 


0 and 


If MxMy^'^ — MyMx^'^ < 0, from f{x 2 ) > 0, /(I) = 0 and f'{l) > 0, f{x) = 0 has a root in [x2,l). 
Similarly, if M^M^^ — MyM^^^ > 0, we have /'(I) < 0, which implies that f{x) = 0 has a root in (1, X3]. 
Also, 1 is not a repeated root of /(x) = 0 since /'(I) / 0 when My / 0. □ 


Proof. ofLemmal4.5l 1. Suppose f{z) = 0. From equation ()4.8p . we have F{z) = 0. So we can 
write F{x) = (x — z)G{z}. Similarly, since Di{z) = 0 (Recall Di{x) = 6^(x) — 4a(x)c(x)), we can write 
Zli(x) = (x — z)D*{x), where D*{x) is a polynomial. It follows that /(x) = (x — z)T{x), where 

T{x) = a{x)[ai{x)f |(x - z)[G{z)f - • 

Since the random walk has genus 1, 2; is not a repeated root of Di{x) = 0, which implies a(x)[ai(x)]^ Aa\x) ^ 
0 (note that a{z) / 0 since Di{z) = 0 and h{z) > 0 when z < 0). It follows that T{z) / 0, that is, 2 : is 
not a repeated root of /(x) = 0. 

2. This is a direct consequence of equation (|4.8p . 

3. Suppose x' is a common root. If a{x') 7^ 0, it is easy to obtain that x' is a branch point. 

Assume a(x') = 0. Clearly, x' cannot be a positive number. Since /i(x') = ai(x')[—26(x')] = 0, 
/o(x') = ^ + 61 (x') = 0 and h{x') / 0, we obtain oi(x') = 0 and 61 (x') = 0, which implies that 

x' = 0 since 61 (x) has only nonnegative zeros. 

4. Let —\z\ be a negative root of /o(x) = 0 in [—X3,—1). From the definition of /o(x), we have 
Y.i>-i,j>oPi^j[-\^\f^di-\z\) = 1, which implies /o(H|) > 0 since yo(|2;|) > |yo(-|2;|)|. According to 
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/o(l) < 0 and Lemma SiH-l, fo{x) = 0 has a root, say z' in (1, \z\). Again, from Yq{\z'\) > |lo(~k^|)|) 
fo{—\z'\) < 0, which implies fo{x) = 0 has a root in {—\z'\,—l) since /o(—1) > 0. Clearly, this root is 
greater than —\z\. 


5. Let |x| G (1,X3]. Since -b{-\x\) < 0, yo(-kl) = = 

f}{—\x\) > ^(|a;|), Di{—\x\) > L>i(|x|) and |c(—|x|)| < c(|x|), we obtain |yo( 


2c(-hl) 


2c(|3:|) 


From 


— \x\ 


< 


K\A)+y/ Di{\x\) 


□ 


Proof, of Lemmal4.6l Assume \z\ = 1 and z 7 ^ 1 or —1. From Lemma 13.21 1. |lo(' 2 )| < 1- Since 


hi{x,y) = ai{x)y+ bi{x) = x \ ^ pijx'y^ - 1 


and when |z| = 1, |i| = 1 as well, we obtain 


*>-id>o 


< X] < 1, 


which yields fo{x) = hi{z,Yo(z}) / 0. 

For z = —1, |yo(—1)1 < 1 if pij is not X-shaped, and yo(—1) = —1 if Pij is X-shaped and p^j is not 
X-shaped. It follows that /o(—1) > 0 in both cases since 6i(—1) > |ai(—1)| and 6i(—1) > 0 in the first 
case and 61 (—1) > |ai(—1)| in the second case. □ 


For special ransom walk 1, we have 


ai(x) = + p^Z\y^ bi{x) = —x 

a{x) = pQ_ix,b{x) = p-ifi — X + pi^qx^ and c{x)=pq^-ix. 
In this case, f{x) becomes 

f{x) = x‘^f*{x), 

where 


(A.l) 

(A.2) 

(A.3) 


f*{x) = dlx^ + dlx-^ + d* 2 X^ + dlx + d, 


= po,ix[l - pZqxZ + [p-1,0 - 3 : + Pi,oxZ [1 - PifiX] + po-ix\pZi + p^ix] 




with 

do = P-i,oPo} and dl = 

Proof, of Proposition UTJ 1. Obviously, From equation (IA.3P and Lemma 14.41 f{x) = 0 has at 
least four real roots with two in [x2,3;3] and two equal to zero. The facts that f{xi) > 0, /(X 4 ) > 0 and 
/(±oo) = —00 yield one root in (—oo,xi] and another root in [x4,+oo). 

2. It is a direct result of Proposition 14. Il l and Lemma 14.51 4. □ 


For the random walk considered in Theorem 14.51 2 (or both pij and p^Z are X-shaped, we have 


ai{x) = pZli + Pljx"^, bi{x) = -X, 


(A.4) 
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a{x) = p-i^i + pi^ix'^, b{x) =-X, c{x)=p-i-i+pi-ix‘^. 

Therefore, f{x) becomes 

f{x) = a{x)bi{x) — b{x)bi{x)ai{x) + c{x)al{x) 

= +Pi,ix^] - 1 + [p-i-i +Pi-ix^]\P-l,i + 

= d^x^ + d^x'^ + d2x‘^ + do, 


(A.5) 


(A. 6 ) 


where 


d% = 


.( 1 ) 

PlA 


1 2 




dA = Pi,i - Pi J + 2pi_ip^}_]^l + p_i-1 


.( 1 ) 

Pi,i 


1 2 


d2 = P-1,1 - Pli,i + 2p_i_ip\[.^p\[ + Pi _1 
1 2 

P-1,-1- 


.( 1 ) 

^- 1.1 


do — 


Ji) 

P- 1,1 


Proof, of Lemmal4.7l /(I) = /(—I) = 0 follows from Tj(l) = —Tj(—1), ai(l) = ai(—1) and 
6i(l) = —6i(—1). From Lemma [4.51 there exists an xq ^ 1, xq £ [x 2 -,xo] such that /(xq) = 0. We 
provide details for the case of fi{z) = 0 and a similar proof can be found for the other case. Since 
X 2 < xo < xo, —b{z) = —b{—z) > 0. It follows that 


-b{xo) , \/b‘^{xo) - 4a(xo)c(xo) 

^ - 

and 

/i(xo) = ai(xo)Pi(xo) + 6i(xo) = 0. 

On the other hand, from —X 3 < —xq < — 1 we have — 6 (—xq) < 0, which yields 


Y(_ ^(^ 0 ) y/b‘^{xo) - 4a{xo)c{xoj 

^ 2a(xo) 2a(xo) 


-’^i(a:o) 


and 

fii-xo) = oi(-xo)yi(-xo) + &i(-xo) = -/i(xo) = 0 . 
It follows from equation (|A. 6 p that /(x) can be written as 


/(x) = do{x‘^ - l)(x^ - Xo)(x^ + p). 


Since ^ > 0, we have p > 0, which indicates that /(x) 


0 has two complex roots. 


(A.7) 

(A. 8 ) 


(A.9) 


□ 


Proof, of Proposition 14.21 Suppose that one of the two complex roots is a root of /o(x) = 0. First 
assume ^ < 1. Then, = ^ < 1 implies |p| < 1. In the case of ^ > 1, we choose a path £ to 
connect the random walk here to the one with ^ < 1. Then on £, the two complex roots of /(x) = 0 
have to pass through the unit circle, which is impossible according to Remark 14.71 and Lemma 14.61 
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